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1  Introduction 


In  this  report  we  present  the  work  performed  in  the  period  May  1990  April 
1992  under  contract  AFOSR-90-022S:  Mathematics:  Numerical  solution  of 
inverse  problems  in  acoustics,  Principal  Investigators:  Luciano  Misici  and 
Francesco  Zirilli. 

In  section  2  we  present  the  work  done  in  the  contract,  in  section  3,4  and  5 
three  scientific  papers  that  explain  in  detail  the  work  discussed  in  section  2 
are  reproduced. 


2.  Statement  of  the  work  accomplished 

Let  R3  be  the  three  dimensional  euclidean  space,  D  C  R3  be  a  bounded 
simply  connected  domain  that  contains  the  origin  with  boundary  dD.  Let 
ii'(x)  be  an  ’’incoming”  plane  wave,  that  is: 

u'(x)  =  e,k(^  (2.1) 


where  k  >  0  is  the  wave  number,  a  is  a  given  unit  vector  and  x  =  (x,  y,  z )  € 
R3.  Let  u3(x)  be  the  scattered  acoustic  field  and  u(x)  be  the  total  field  that 
is: 

u(x)  =  u‘(x)  +  u3(x )  (2-2) 

Let  us  consider  the  following  ’’direct”  problems: 

Problem  2. 1-  Find  u  defined  in  R3  \  D  such  that 


\0u 3 
dm  r  —— 

r - OO  l  (JT 


0  in  R 3  \D 

(2.3) 

on  dD 

(2.4) 

—  iku3 1  =  0 

(2.5) 

where  A  =  and  r  =  { x2  +  y2  +  z2)lf2. 


Problem  2.2-  Find  u  defined  in  R3  \  D  such  that 


A u  4-  k2u  =  0  in  R 3  \  D 


(2.6) 


on  dD 


r  du3 

hm  r  — — 
r— oo  or 


(2.7) 

(2.S) 
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where  u  is  the  outward  normal  vector  to  dD. 

(2.9) 

(2.10) 

(2.11) 

Problem  2.1  and  in  particular  the  Dirichlet  boundary  condition  (2.4) 
corresponds  to  the  study  of  acoustically  soft  obstacles,  Problem  2.2  and 
in  particular  the  Neumann  boundary  condition  (2.7)  corresponds  to  the 
study  of  acoustically  hard  obstacles,  finally  Problem  2.3  and  in  particular 
the  mixed  boundary  condition  (2.10)  correspomds  to  the  study  of  obstacles 
characterized  by  an  acoustic  impedance.  Moreover  we  have: 

u3(x)  =  —  F(x,k,a)  +  0(\)  ;  r  — >  oo  (2.12) 

r  ri 

where  x  =  ^  and  F(i,k,a)  is  the  far  field  pattern  associated  to  u3(x) 
generated  by  u‘(x).  Let  An,  n  =  1,2. ...  be  the  eigenvalues  of  the  Helmholtz 
equation  (2.3)  considered  in  D  with  the  appropriate  boundary  conditions. 
Let  B  =  {x  €  R3  |  ||x||  <1}  and  dB  =  {x  6  R3  |  ||x||  =  1},  we  will  consider 
the  following  inverse  problems: 

Problem  2.4-  Given  k  >  0,  —  k2  ^  A„,  n  =  1,2,...  and  Q  dB  from 

the  knowledge  of  the  far  field  patterns  F(x,k,a )  of  problem  2.1  (Dirichlet 
boundary  value  problem)  for  a  €  12 1  ,  x  G  Q2  find  the  boundary  of  D,  dD. 

Problem  2.5-  Is  analogous  to  Problem  2.4  when  Problem  2.2  that  is  the 
Neumann  boundary  condition  is  considered  instead  of  Problem  2.1. 


Problem  2.3-  Find  u  defined  in  R 3  \  D  such  that 


Au  +  k2u  =  0  in  R 3  \  D 


du 

lliU  +  a o-r-  =  0  on  dD 
au 


lim  r 
r— ■  00 


du3 

dr 


—  iku3 


=  0 


where  p.\  and  pi  are  given  constants. 
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Problem  2.6-  Is  analogous  to  Problem  2.4  when  Problem  2.3  that  is  the 
mixed  boudary  condition  is  considered  instead  of  Problem  2.1. 

The  inverse  Problems  2. 4, 2. 5, 2. 6  have  been  considered  with  special  attention 
to  the  resonance  region  case,  that  is  the  case  when 

kL  ~  1 

where  L  is  a  characteristic  lenght  of  the  obstacle. 

In  the  research  reported  here  we  have  considered  two  types  of  obstacles: 

(1)  obstacles  with  smooth  boundary  dD.  For  this  class  of  obstacles  have 
developed  a  numerical  method  to  solve  the  inverse  Problems  2. 4, 2. 5, 2. 6 
(see  the  references  reported  in  sections  3,4,5  of  this  report).  The 
numerical  method  developed  generalizes  the  method  introduced  by 
Colton  and  Monk  and  uses  the  Herglotz  function  technique. 

(2)  obstacles  with  Lipschitz  continuous  boundary.  For  this  class  of  obstacles 
we  have  studied  [1]  the  direct  Problems  2. 1,2. 2, 2. 3.  The  existence  and 
uniqueness  of  the  solution  of  problems  2. 1,2. 2, 3. 3  has  been  established. 
In  order  to  compute  the  far  field  patterns  generated  by  these  non 
smooth  boudaries  we  have  generalized  the  method  suggested  by  Milder 
in  [2].  However  up  to  now,  due  to  the  difficulties  encountered  in 
computing  the  far  field  data,  we  have  not  been  able  to  work  out  a 
satisfactory  method  to  solve  the  inverse  Problems  2. 4, 2. 5, 2. 6  when  non 
smooth  boundaries  are  considered. 

References 

[1]  A.  Vecchio,  F.  Zirilli,  in  preparation. 

[2]  D.M.  Milder:  An  improved  formalism  for  wave  scattering  from  rough 
surfaces,  J.  Acoust.  Soc.  Am.,  89,  1991,  529-541. 
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L.  Misici,  F.  Zirilli:  ’’An  inverse  problem  for  the  three  dimensional 
Helmholtz  equation  with  Neumann  or  mixed  boundary  conditions:  a 
numerical  method”,  in  Mathematical  and  numerical  aspects  of  wave 
propagation  phenomena,  G.  Cohen,  L.  Halpern,  P.  Joly  Editors,  SIAM 
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An  Inverse  Problem  for  the  Three  Dimensional 
Helmholtz  Equation  with  Neumann  or  Mixed 
Boundary  Conditions:  A  Numerical  Method’1' 

Luciano  Misicif 
Francesco  Zlrilli t 


Abstract.  A  numerical  algorithm  for  a  three  dimensional  inverse  acoustic  scattering 
problem  is  considered.  From  the  knowledge  of  several  far  field  patterns  of  the  Helmholtz 
equation  a  closed  surface  3D  representing  the  boundary  of  an  unknown  obstacle  D  is 
reconstructed.  The  obstacle  D  is  supposed  to  be  hard  (i.e.  Neumann  boundary  condition) 
or  characterized  by  an  acoustic  impedance  (i.e.  mixed  boundary  condition). 

1.  Introduction 

Let  Ra  be  the  three  dimensional  euclidean  space,  £  =  (r,y,r)  6  HP  be  a  generic  vector, 
(...)  will  denote  the  euclidean  scalar  product  and  ||.||  the  euclidean  norm.  Let  D  C  R3 
be  a  bounded  simply  connected  domain  with  smooth  boundary  dD  that  contains  the 
origin.  Let  u‘ (z)  be  an  incoming  acoustic  plane  wave,  that  is: 

u'(x)  =  «•'*<**)  (1.1) 

where  k  >  0  is  the  wave  number  and  a  €  R3  is  a  fixed  unit  vector.  Let  us  denote  with 
u‘(z)  the  acoustic  held  scattered  by  the  obstacle  D  and  with  u(z)  the  total  acoustic  field, 
that  is: 

«(S)  =  «•'(!) +  «*(S)  (1-2) 

The  total  acoustic  field  u(z)  satisfies  the  Helmholtz  equation: 

:«ri!»  =  0  in  R 3\D  (1.3) 

and  the  scattered  acoustic  field  u’(z)  satisfies  the  Sommerfeld  radiation  condition  at 
infinity,  that  is: 

du* 

lim  r (  — - ikn'}  =  0  ( 1 . -t ) 

r  —  co  dr 

The  research  reported  in  this  paper  Iia.i  been  made  possible  through  the  support  and  sponsorship 
of  the  Uaited  States  Government  through  the  Air  Force  Office  of  Scientific  Research  under  contract  n^ 
AFOSR  90-0228  with  the  University  di  Cnmerino 

f Dlpartlmento  dl  Matematlc*  e  Fisica  University  di  Cameriuo  >82032  C&merino  (UC)-  Italy 
fDIp  artiraento  dl  Xdatematiea  *  G.  Cajtelnuovo*  University  dl  Roma  "La  SApienta"  >00186  Romv 

ItaJy 


3 


49S 


Misici  and  Ziriili 


where  A  =  j~:  +  +  J~r  <»  the  laplacian  and  r  =  i|r|j  =  (r3  +  y3  -r  r3)1 1,3  .  Moreover  the 

total  acoustic  Geld  u(r)  satisfies  a  boundary  condition  on  dD.  This  boundary  condition 
can  be  formulated  in  several  different  ways,  depending  on  the  nature  of  the  obstacle  D. 
In  (l],[2|  we  have  considered  the  acoustically  soft  obstacles  [3|  that  are  characterized  by 
the  Dirichlet  boundary  condition: 

u  =  0  on  d  D  (1.5) 

In  this  paper  we  restrict  our  attention  to  the  acoustically  hard  obstacles  [3]  characterized 
by  the  Neumann  boundary  condition 


—  =  0  on  dD  (1.6) 

a  v 

where  if  is  the  unit  normal  on  dD,  and  to  the  obstacles  characterized  by  an  acoustic 
impedance  x  [3|  that  satisfy  the  mixed  boundary  condition 

du 

u  +  y—  =0  on  dD  (1.7) 

dir 

We  assume  that  x  is  a  given  constant.  We  consider  two  boundary  value  problems:  the 
Neumann  boundary  value  problem  given  by  (1.3),(I.4),(1.6)  and  the  mixed  boundary 
value  problem  given  by  (1.3),(1.4),(1.7).  In  [4|  it  is  shown  that  the  scattered  Geld  u’(z) 
of  the  Neumann  and  mixed  boundary  value  problem  has  the  following  expansion 


u'(x)  =  F0(=,k,a)  +0(  — )  when  r  -  co 


where  i  =  jjfjj  ,  a  5=  0  and  Fn{i,k, a)  is  the  far  Geld  pattern  generated  by  the  incoming 
wave  (1.1)  that  hits  the  obstacle  D. 

In  this  paper  we  introduce  a  numerical  method  for  am  inverse  problem  for  the  three 
dimensional  Helmholtz  equation,  that  is  from  the  knowledge  of  the  nature  of  the  obstacle 
i.e.  the  boundary  condition  on  dD  satisfied  by  u(;)  ana  from  the  far  Gelds  Fa  generated 
by  several  incoming  waves  we  want  to  recover  the  shape  of  the  obstacle  dD. 

The  inverse  acoustic  scattering  problem  have  received  a  lot  of  attention  in  the  scientiGc 
and  technical  literature,  here  we  wiil  refer  only  to  the  work  of  Colton  and  Monk  [5|,[6),|7] 
since  the  work  presented  in  this  paper  has  beeen  inspired  by  (7).  To  be  more  precise 
let  An  ,  n  =  1,2,...  be  the  eigenvalues  of  the  Hcimhoitz  equation  in  the  interior  of  D, 
with  Neumann  boundary  condition  { 1  5)  or  with  the  mixed  boundary  condition  (1.7);  let 
B  =  (x  €  R3  |  Usd  <  1),  and  33  be  the  boundary  of  3.  We  will  consider  the  following 
inverse  problem: 

Problem  1 . 1  Let  us  assume  that  *u'r)  satisfies  the  Neumann  boundary  condition  (1.6)  or 
the  mixed  boundary  condition  ( 1.7).  Let  TZ ,  C  d B  ,  Hj  C  {£  €  R  |  <f  >  0}  be  two  given 
sets  such  that  A,  €  D-i  •  =  1,  2...  From  the  knowledge  of  Fa{£ ,  k,a),  for  a  €  Oi  ,  k  € 
determine  the  boundary  of  the  obstacle  dD. 

We  note  that  the  condition  A,  fl-j  f  =  1,2...  is  a  non-resonance  condition,  that  Hi  is 
the  set  of  the  directions  of  the  incoming  waves  and  that  the  far  field  F„  is  observed  for 

Z.  €  d  B. 

In  this  paper  we  present  a  numerical  method  to  solve  Problem  1.1,  in  particular  in  section 
2  we  derive  the  mathematical  relations  needed  to  develop  our  method,  in  section  3  we 
present  our  numerical  method,  finally  in  section  4  some  numerical  experience  is  shown. 
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2.  The  mathematical  formulation  of  the  inverse  problem 
For  €  R3  let 


*(*lls-yll)  = 


e**ll*-yll 

Its  -  yj| 


be  the  Green’s  function  of  the  Helmholtz  operator  with  the  Sommerfeld  radiation  condi¬ 
tion  at  infinity.  It  is  easy  to  see  that: 

*Wls-S£lI)  =  ^«-,'*,^,+0(^j5)  when  ||=i|-=o  (2.2) 

moreover  from  the  Helmholtz  formula  [8]  we  have: 

where  da (y)  is  the  surface  measure  on  3D. 

Substituting  (2.2)  in  (2.3)  and  using  (1.8)  we  have: 

Fa[Lk,a)  =  ~(  (u(y)^,  -  (2.4) 

4*  Jao  ~  at/(y)  du(g) 

Let  g(i )  6  L^(dfl,  dc)  where  L '(3  B ,  da)  is  the  space  of  square  integrable  functions  with 
respect  to  the  measure  da  and  y(jj)  is  the  complex  conjugate  of  j(i). 

For  every  g{i)  g  L 2{3B,da)  from  (2.4)  interchanging  the  integrals  we  have: 

Intern 

1  f  dv(Jcy)  3u(y) - 

‘zLw-sbg 


waere 

«(*v)  =  [  9{iVk[kv-]da{k)  (2.6) 

Job 

It  is  easy  to  see  that  v(ky)  satisfies  the  Helmholtz  equation  for  y  g  R3,  moreover  v(ifcy) 
is  the  Herglotz  wave  function  corresponding  to  the  Herglotz  kernel  y(r).  Since  the  total 
acoustic  field  u  satisfies  the  boundary  condition  (1.6)  or  (1.7)  on  the  surface  3D  formula 
(2.5)  reduces  to: 

(•)  /  =  7"  f  u(y)TTTT<^(y)  (2-7) 

JOB  Jod  ~  °V\3L) 

when  the  Neumann  boundary  condition  (1.6)  is  satisfied,  or  to 

(.'o  =  ± (2.3) 

when  the  mixed  boundary  condition  (1.7)  is  satisfied. 

We  restrict  our  attention  to  the  Neumann  Herglotz  domains,  that  is  domains  such  that 
f’  unique  solution  u  of: 

(A  +  Jc2)u  =  0  z€D  (2.9) 


! 
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,-MUlI 


z<=3D 


du  d»{z)  JfeiUH  -  ~  v  '  7 

is  given  by  (2.6)  for  a  suitable  choice  gn{z,  k)  of  g(z),  or  to  the  mixed  Herglotz  domains, 
that  is  domains  such  that  the  unique  solution  v  of 

(^  +  k2)v  =  0  ;S  D  (2.11) 

d  d 

^  +  +  (2-12) 
is  given  by  (2.6)  for  a  suitable  choice  gn[z,  k)  of  g(x). 

We  note  that  in  the  definition  of  Hergloti  domains  we  have  exploited  the  hypotesis 
k'2  jt  X,-  i  =  1,2,  •••  .  A  simple  computation  shows  that  the  sphere  is  an  Hergloti 
domain  that  is  the  claas  of  the  Hergloti  domaina  ia  not  empty.  In  (2. 7), (2. 8)  let  v  be  the 
Herglotr  wave  function  associated  to  gu(d,k),  using  (2. 10), (2. 12)  respectively  and  the 
Helmholtz  formula  (2.3)  we  have 


f  ?„{=,*,&)  9  hO.,*)  do  (z)  = 

J  BB  * 


when  Jw(Si^)  iJ  the  Herglotz  kernel,  formula  (2.13)  holds  Vfc,  q.  Problem  1.1  proposed 
in  section  1  will  be  solved  in  three  steps: 

(i)  from  th.  knowledge  of  some  far  fields  F„  using  (2.13)  determine  the  Herglotz  kernel 
?tf(ii£)  of  the  domain  D 

(ii)  from  gff{z,k )  using  (2.6)  find  the  corresponding  Herglotz  wave  function  v 

(iii)  determine  3D  using  (2.10)  (Neumann  boundary  condition)  or  (2.12)  (mixed  boun¬ 
dary  condition). 

3.  The  numerical  method 

Given  D  C  R3  and  the  boundary  conditions  (1.5)  or  (1.7)  satisfied  by  u  on  3D,  let 
0,  =  (o,  €  dB  |  f  =  1 , 2 ,  •  -  * ,  >V }  be  the  set  of  directions  of  the  incoming  waves  and 
flj  =  {£;■  >  0  |  f  =  1, 2,  •  •  • ,  M}  be  the  set  of  wave  numbers  of  the  incoming  waves.  The 
data  of  our  problem  are  the  measurements  of  the  far  fields  F0(i,  k,  ,au)  i  £  d B  ;  g  — 
1, 2,  ••  • ,  M  ;  i  =  l,2 ,  •  •  • ,  W.  In  the  numerical  experience  shown  in  section  4  the  data  are 
obtained  by  solving  numerically  the  "direct”  problems  (1.3), (1.4), (1.6)  or  (1.3), (1.4), (1.7). 
Let  (9,<f>)  be  the  polar  angles  so  that 

z(9 ,  <?)  =  (sin  9  cos  <p,  sin  9  sin  d>,  cos  9)  (3-1) 


f/im(z)  =  7lm/5i'n(co3  9)  COS  m<p  ,  V'fm(i)  =  7/mP/m(cos  9)  sin  m4> 

be  the  spherical  harmonics  that  is  Pi  are  the  Legendre  polinomials,  Pt"'  is  the  mtK 
derivative  of  Pi  and  nm  are  the  normalization  factors  in  L7{3B,do).  Prom  these  data 
our  computation  proceeds  in  four  steps: 

Step  1.  For  each  j  compute  the  Fourier  coefficients  of  F0(s.,  k},^J  ,  * ’  —  1,2,-  •  - ,  iV . 
Given  >  0  we  assume  that  the  far  field  F0(z,  k}  ,a,)  can  be  approximated  by  a 

truncated  Fourier  series,  that  is: 

F.&ki.Oi)  =  E  +  E  E  F'^vUi)  (3-2) 
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The  Fourier  coefficients  {F^lt  Ft'^2}  are  determined  computing  the  following  integrals: 


r;L  =  J 

f  F0{z,k.;&)Ulm(z)d*(z) 

33 

,  i  =  0,1, 

•«»  ^rrxax 

;  m  =  0,  1, ..,/ 

(3.3) 

f;L  =  j 

f  F0(z,k]l&A)Vlm(z)da(z) 

OB 

,  /=1.2, 

■  •  i  ^m<u 

;  m  =  1,2,  ..,/ 

(3.4) 

Step  2.  From  the  coefficients  of  F0[z,k;  ,a,)  ;  i  =  1,2,  ...,1V  to  the  Herglotz  kernels 
gn{z,kj)  ,;=  1,2 of  the  domain  D. 

Let  gn(s.<ki)  £  L'(33,da)  be  the  Herglot:  kernel  associated  to  the  domain  D  when  we 
use  the  wave  number  kj  £  n2.  We  assume  for  gn (i,  k..)  the  expression  of  a  truncated 
Fourier  expansion,  that  is 


i*. 


y*  U.  *y)  =  EE  sL  M=)  +  E  E 


(3.5) 


is=0  m  =  0 


f=s  1  m=s  1 


where  0  <  L7  <  From  (2.13)  using  the  orthogonality  properties  of  the  spherical 

harmonics  we  have: 

_ _  £  _ 

d*  E]  E  ^trrxl^im.2  ~ 

l  =0  m  =  0  J  =  1  m-—  1 


7?-  ,  •  =  1,2 . N  ;  ;  =  l,2t...tAf  (3.6) 

k i 


So  that  for  each  j  the  Fourier  coefficients  are  determined  solving  (3.6).  For 

j  fixed  the  linear  system  (3.6)  has  IMA* t?.l  4 •  unknowns,  that  is  in  order  to 

determine  we  need  N  incoming  waves  with  iV  >  +  4- 

If  the  obstacle  D  has  some  symmetry  such  as  cylindrical  symmetry  around  the  z-axis 
and/or  symmetry  with  respect  to  the  equator  then  a  similar  symmetry  can  be  assumed 
on  gn{k ,ki)-  This  assumption  reduces  substantially  the  number  of  unknowns  in  (3.6) 
and  as  a  consequence  the  number  N  of  incoming  waves  needed  to  recover  the  desired 
approximation  of  gn{k >kj)  (2). 

Step  3.  From  the  Herglotz  kernels  gn{k>’K])  to  the  Herglotz  wave  functions  v;(l:Jy). 
From  (2.6)  and  (3.5)  we  have: 


v}(k}y)  =  4*  (E  E  s/mi1' ^(MyllMmfy)  +  E  E 


A(M*ll)VUi)  }  (3-7) 


i  JsO  m  =  0 


Jr:  l  rnss  1 


where  y  =  jj^jj  and  J\(-)  is  the  spherical  Bessel  function  of  order  l. 

Step  4.  From  the  Herglotz  wave  functions  vj[k}y)  ;  j  =  1,2,...,M  to  the  boundary  of 
the  obstacle  3D- 

Let  [r,9,$)  be  the  polar  variables  we  assume  that  exist  0  <  a  <  &  <  oo  and  a  function 
/(M)  with  a  <  /  <  6  such  that  3D  =  (r  =  f[8,$)  |0<5<a-,0<<^<  2  x).  We 
approximate  f(9,<l>)  with  a  truncated  Fourier  series,  that  is: 

L.  i  L,  ( 

/{M)  =  E  E  C/-^/m(i)+EE  (3-8) 

1=0  m  =  0  1=1  m=  1 

where  L,  >  0  is  chosen  depending  on  3D  i.e.  simple  obstacles  can  be  reconstructed  with 
Lf  =  4,6.  Moreover  if  3D  has  some  symmetries  these  can  be  trnslated  in  properties 
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of  the  coefficients  {c/m i,c(rnl}.  Let  c  =  {ciml,ctm3}  0  <  l  <  L„  ,  0  <  m  <  /  be  an 
— ILL*— d  +  LA ^ t 1  i  dimensional  vector.  The  unknown  boundary  dD  is  obtained 
minimizing  with  respect  to  c 


dl,  1  f J,r  /"  I  3  e~'k’f  2 

**L  +  3ia9d 9  (3-9) 

when  the  Neumann  boundary  condition  (1.6)  is  considered,  or 

=  d*fo  \(xjz  +  lKu’  +  Lkj-]\ 3ia9d9  (3-10) 

when  the  mixed  boundary  condition  (1.7)  is  considered. 

In  (3.9), (3. 10)  the  Herglotz  wave  function  u;  are  computed  In  (r  =  /(0,^),0,d)  and  v  is 
the  unit  exterior  normal  to  the  surface  r  =  J{9,4>)  and  /  is  given  by  (3.8).  The  integrals 
in  (3.9), (3. 10)  are  approximated  with  some  elementary  quadrature  formula.  The  weights 
are  introduced  in  order  to  make  the  different  terms  in  the  sums  over  j  of  the  same 
order  of  magnitude.  When  the  minimization  of  the  functions  I,  (c)  ,  i  =  1,2  does  not 
give  a  satisfactory  reconstruction  of  dD  we  minimize 

P.(c)  =  /,(<:)  +  u,(c)  ;  ,-=1.2  (3.11) 

where  w(c)  is  a  penalization  term. 

For  a  large  class  of  surfaces  <5  D  including  the  ones  with  cylindrical  symmetry  with  respect 
to  the  z-axis  we  have  |^  =  ££=Oat0  =  O  (North  pole)  and  0  =  x  (South  pole),  moreover 

if  dD  is  also  symmetric  with  respect  to  the  equator  we  have  =  0  at  8  ~  (the 

Equator). 

When  =  0  the  relations 


(3.12) 


and 

3  e~'k>f 

(*  ^  +  i)K  + -j-y)  =  0  (3.13) 

become  nonlinear  equations  in  the  unknown  /  that  can  be  solved,  for  example,  using 
the  bisection  method.  In  this  way  we  can  obtain  ,f3,f3  estimates  of  f{9,4>)  when 
8  =  01,02,9a  with  0j  =0,02  =  The  penalization  term  w(z)  is  given  by: 


3 

w(c)  =  ^Pi(/(0,,d)-/;)J  (3.14) 

1=  I 

where  /(0,d)  ■*  given  by  (3.8)  and  p;  >  0  ,  i’  =  1,2,3  are  weight  factors.  The  mini¬ 
mization  of  the  functions  /,(c)  or  P,(c )  ,  «  =  1,2  is  performed  with  a  stochastic  global 
minimization  algorithm  introduced  in  (9j,(i0j. 


4.  The  numerical  experience 

In  our  numerical  experience  we  chosen  in  (1.7)  the  acoustic  impedance  x  —  1-  The 
surfaces  3D  considered  are  the  following  ones: 
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1- 

Oblate  Ellipsoid 

(§*)2+(2-y)2  +  x2  =  l 
x3  +  y2  +  (fr)2  =  1 

(4.1) 

2- 

Prolate  Ellipsoid 

(4.2) 

3- 

Short  Cylinder 

((!*)3  +  (§yn6+*l0  =  i 
(**  +  » 3)s  +  (!^),0  =  i 

(4.3) 

4- 

Long  Cylinder 

(4.4) 

5- 

Vogel’s  Peanut 

r  =  |  (cos2  5  +  4  sin2  0)1/2 
r  =  jf(Jf  +  2C0S35)1/2 

(4.5) 

6- 

Pseudo  Apollo 

(4.6) 

All  these  surfaces  axe  cylindrically  symmetric  with  respect  to  the  z-axis  and  the  surfaces 
1,2, 3,4, 5  are  also  symmetric  with  respect  to  the  equator.  These  symmetries  are  always 
exploited  to  reduce  the  number  of  Fourier  coefficients  in  the  expansions  of  the  Herglotz 
kernels  and  of  the  surfaces  f[9,4>). 

We  observe  that  the  obstacles  D  corresponding  to  1,2, 3, 4  are  convex  and  the  ones  corre¬ 
sponding  to  5,6  are  not  convex.  Finally  a  characteristic  length  L  of  the  obstacles  can  be 
chosen  equal  1.  The  set  of  the  directions  of  the  incoming  waves  is: 

ni  =  {(«;. o)  i  a;  ,  /  =  o,  i . if)  (4.7) 

with 

N  =  Lmax  +  1  (4.8) 

The  set  fls  is  a  subset  of  the  set  (2,3)  and  Lg  =  We  observe  that  with  this  choice 

of  Q-2  the  product  kjL  is  of  order  one,  that  is  we  are  working  in  the  resonance  region. 
For  j  =  0, 1,  •  •  • ,  36  let  9j  =  ^  ,  /(0y,O)  be  the  exact  value  of  the  surface  given  by 

(4.1) . (4.6)  and  fc(9 y,0)  be  the  value  reconstructed  performing  the  numerical  procedure 

described  in  section  3.  The  relative  L2  error  in  the  points  {$;  j  j  =  0, 1,... ,36),  that  is 


Ci:^o(/(g|.0)-/.(g;.0))a 

I  r-'-JS  tlln  n) 


1/2 


(4.9) 


is  used  as  a  performance  index. 

The  results  obtained  are  shown  in  Table  4.1,  Table  4.2,  Fig.  4.1a,  4.1b,  4.1c,  Fig.  4.2a, 
4.2b,  4.2c  and  Fig.  4.3a,  4.3b,  4.3c. 


Table  4.1  Hermann  problea 


Object 

l 

Recons  tructicrJL^aj 

*2 

Penaiizati 
North  Pole 

on  term 
Equator 

Oblate  Ellipsoid 

1  1  3  |  4 

2 

nc 

yes 

yes 

p.00415 

m 

2  3  14 

213 

no 

no 

D  .004 13 

Prolate  Ellipsoid 

1 

yes 

yes 

MMiiMi 

n 

4  3  1  4  1  2  |  3 

no 

no 

.Iil.kHI 

isM  J  -  I*. 

Short  Cylinder 

5  ■  8  |  5  i  2  Snc 

yes 

yes 

0.0476 

n 

6  |  3  i  5 

2  f  3 

no 

no 

0.0361 

Long  Cylinder 

7  i  3 

6 

2  ha 

yes 

yes 

0.0551 

t • 

3 

3 

6 

2 

no 

□  0 

0.0785 

Vogel’s  Peanut 

9 

8 

< 

2 

nc 

yes 

yes 

0.0299 

n 

10 

3 

4 

3J 

no 

no 

0.0225 

Pseudo  Apollo 

11 

8 

Q 

B 

£ 

yes 

failure 

ft 

12 

8 

a 

B 

B 

no 

no 

failure 

ft 

13 

8 

a 

B 

B 

yes 

no 

0.0353 

14 


504 


MLsici  and  Zirilli 


Table  4.2  Problaa  with  alxad  boundary  condition 


Object 

■ 

H 

I 

1 

WBMIDL 

Oblate  Ellipsoid 

1 

8 

1 

E 

no 

yes 

yes 

0.0102 

If 

... 

2 

8 

□ 

B 

B 

no 

no 

3 

8 

< 

2 

yes 

yes 

0.0139 

m 

4 

8 

□ 

B 

B 

no 

no 

0.0232 

Short  Cylinder 

5 

s 

B 

B 

£ 

yes 

yes 

0.0399 

n 

6 

8 

B 

B 

B 

no 

no 

0.0568 

Long  Cylinder 

7 

8 

B 

B 

£ 

yes 

yes 

0.0424 

n 

8 

8 

B 

B 

B 

no 

no 

0.0293 

Vogel's  Peanut 

9 

8 

□ 

B 

£ 

yes 

yes 

0.0448 

If 

10 

8 

□ 

B 

B 

no 

no 

0.0495 

Pseudo  Apollo 

11 

8 

□ 

B 

£ 

yes 

yes 

failure 

ft 

12 

8 

□ 

B 

B 

no 

no 

failure 

n 

13 

8 

□ 

B 

B 

yes 

no 

0.0356 

Fig.  4.1a  Original 
Long  Cylinder 


Fig.  4.1b  Reconstruction 
n°  8  of  table  4.1 


Fig  4.1c  Reconstruction 
n"  8  of  table  4.2 


Fig.  4.2a  Original 
Vogel's  Peanut 


Fig.  4.2b  Reconstruction 
n”  10  of  table  4.1 


Fig  4.2c  Reconstruction 
n®  10  of  table  4.2 
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Fig.  4.3a  Original 
Pseudo  Apollo 


Fig.  4.3b  Reconstruction 
n°  13  of  table  4.1 


Fig  4.3c  Reconstruction 
n°  13  of  table  4.2 


The  reconstruction  procedure  adopted  here  based  on  the  minimization  of  the  functions 
/i  or  P{  ,  i  =  1,2  is  more  robust  than  the  reconstruction  procedure  adopted  in  |l),[2| 
that  integrates  numerically  an  initial  value  problem  for  an  ordinary  differential  equation. 
However  the  global  minimization  of  the  functions  /;  or  F,  ,  i ’  =  1,2  using  a  stochastic 
algorithm  is  computationally  much  more  expansive  than  the  solution  of  an  initial  value 
problem.  The  reconstruction  technique  suggested  in  this  paper  can  be  adapted  to  the 
problem  with  Dirichlet  boundary  condition  with  very  satisfactory  results. 
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1-  Introduction 


Let  R3  be  the  three  dimensional  euclidean  space,  x  =  (i ,y,z)T  6  R3  be  a  generic  vector  and  the 
superscript  T  denotes  the  transpose  operation.  The  euclidean  scalar  product  will  be  denoted  with 
(•,  •)  and  ||  •  ||  will  denote  the  euclidean  norm. 

Let  D  C  R3  be  a  bounded  simply  connected  domain  with  smooth  boundary  dD,  in  the  following, 
without  loss  of  generality,  we  assume  that  D  contains  the  origin.  Let  u'(x)  be  an  incoming  acoustic 
plane  wave,  that  is: 

u'(x)  =  eifc<£*2>  (1.1) 

where  k  >  0  is  the  wave  number  and  a  6  R3  is  a  fixed  unit  vector  (i.e.  ||a||  =  1).  Let  us  denote 
with  u3(x)  the  acoustic  field  scattered  by  the  obstacle  D  and  with  u(x)  the  total  acoustic  field, 
that  is: 

u(x)  =  u‘(x)  +  u*(x)  (1.2) 

The  total  acoustic  field  u(x)  satisfies  the  Helmholtz  equation: 

Au(x)  +  Aru(x)  =  0  x  e  R3  \  D  (1.3) 

where  A  =  d2 /dx2  +  d2 /dy2  +  d2 /dz2  is  the  Laplace  operator.  The  scattered  field  ua(x)  satisfies 
the  Sommerfeld  radiation  condition  at  infinity 

du3 

lim  r{— - iku3}  =  0  (1.4) 

r— oo  Qt 

where  r  =  |[x||.  Moreover  the  total  acoustic  field  u(x)  satisfies  a  boundary  condition  on  dD.  This 
boundary  condition  is  the  mathematical  counterpart  of  the  physical  character  of  the  obstacle,  that 
is:  for  acoustically  soft  obstacles  we  require  the  Dirichlet  boundary  condition: 

u(x)  =  0  x  6  dD  (1.5) 

for  acoustically  hard  obstacles  we  require  the  Neumann  boundary  condition: 

^>=0  ztBD  (1.6) 

where  n(x),  x  €  dD  is  the  unit  outward  normal  to  dD  at  the  point  x,  and  finally  for  obstacles 
characterized  by  an  acoustic  impedance  we  require  the  mixed  boundary  condition: 

u(x)  +  =  0  xedD  (1.7) 
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We  assume  x  to  be  a  real  constant  to  derive  the  relations  of  section  3,  however  this  assumption 
can  be  avoided. 

We  call  direct  problem  the  problem  of  determining  the  scattered  field  ua(x),  x  €  R3  \  D  given 
the  incoming  field  u'( x)  the  obstacle  D  and  its  physical  character,  that  is  given  one  of  the  three 
boundary  conditions  (1.5), (1.6), (1.7).  The  direct  problem  is  a  boundary  value  problem  for  the 
Helmholtz  equation,  and  has  been  widely  studied.  For  an  exposition  of  several  mathematical  results 
on  the  direct  problem  see  for  example  [1].  In  particular  it  can  be  shown  [1]  that  the  scattered  field 
u3(x)  corresponding  to  the  boundary  condition  (1.5),  (1.6)  or  (1.7)  has  the  following  expansion: 

ua(x)  =  - —  F0(x,k,a)  +  0(-^)  r -*  oo  (1.8) 

so  that  when  r  — *  oo  the  leading  term  of  the  expansion  in  inverse  powers  of  r  is  given  by  a 
spherical  wave  elkr/r  coming  out  from  the  origin  modulated  by  the  ”far  field”  F<j.  We  note  that  F0 
depends  on  k,a  that  are  the  parameters  characterizing  the  incoming  wave  (1.1)  and  on  f  =  -jgjj, 
for  x  /  0.  We  note  that  the  Helmholtz  equation  (1.3)  is  obtained  from  the  wave  equation  assuming 
that  the  incoming  field  and  the  corresponding  scattered  field  are  time  harmonic,  that  is  their  time 
dependence  is  given  by  a  factor  elwt  where  u  is  a  constant. 

The  inverse  problem  that  we  consider  here  is  the  following:  given  the  character  of  the  obstacle  (i.e. 
acoustically  soft,  hard,  or  characterized  by  an  acoustic  impedance)  and  the  far  field  Fo(x,  k,a)  for 
one  or  several  incoming  waves  u‘(x)  with  different  incident  directions  a  and/or  wave  number  k 
determine  the  boundary  of  the  obstacle  dD.  This  inverse  problem  is  known  to  be  ill  posed  and  due 
to  its  great  interest  both  in  mathematics  and  in  several  application  fields  has  been  widely  studied, 
for  a  review  see  [2]. 

The  numerical  methods  used  to  solve  the  inverse  problem  considered  here  can  be  divided  in  two 
types:  the  first  type  consists  of  an  iterative  procedure  that  at  each  step  requires  the  numerical 
solution  of  a  direct  problem,  the  second  type  consists  of  genuine  methods  for  the  inverse  problem 
that  do  not  require  the  solution  of  the  direct  problem.  In  the  first  type  we  mention  the  work  of 
Roger  [3],  Murch,  Tan  and  Wail  [4],  Wang  and  Chen  [5],  Angell,  Colton  and  Kirsch  [6],  Kristensson 
and  Vogel  [7].  In  the  second  type  we  mention  the  work  of  Kirsch  and  Kress  [8],  [9],  [10]  and  the 
work  of  Colton  and  Monk  [11]. 

In  this  paper  we  introduce  a  numerical  method  to  solve  this  inverse  problem  based  on  the  Herglotz 
wave  function  method  introduced  by  Colton  and  Monk  in  [11]  and  further  developed  by  the  authors 
in  [12],  [13],  [14],  [15].  In  particular,  based  on  previous  work  by  the  authors  [14],  we  extend  the 
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Herglotz  function  method  introduced  in  [11]  for  acoustically  soft  obstacles  to  hard  obstacles  or 
obstacles  characterized  by  an  acoustic  impedance.  The  analytical  relations  obtained  are  exploited 
to  built  up  numerical  algorithms.  Finally  these  algorithms  are  efficient  in  the  so  called  resonance 
region,  that  is  when 

kL  =  0(1)  (1.9) 

where  L  is  a  characteristic  length  of  the  obstacle  D.  In  section  2  we  derive  the  analytical  relations 
that  are  the  basis  of  the  numerical  methods.  In  section  3  the  basic  numerical  method  developed 
for  the  solution  of  the  inverse  problem  is  presented.  In  section  4  some  special  features  of  the 
reconstruction  procedure  in  the  case  of  acoustically  soft  obstacles  are  shown.  Finally  in  section 
5  we  introduce  the  test  problems  used  to  test  the  methods  of  sections  3  and  4  and  we  show  the 
numerical  results  obtained. 


2-  The  mathematical  formulation  of  the  inverse  problem 


For  x,y  6  R3  let 


*(s,y)  = 


-  4*||*  -  *|| 

be  the  Green’s  function  of  the  Helmholtz  operator  that  satisfies  the  Sommerfeld  radiation  condition 
at  infinity.  It  is  easy  to  see  that: 

piiilril  1 

=  +  °(iiiF*  when*-“  <2-2> 

moreover  from  the  Helmholtz  formula  [14]  we  have: 

f  r  ,  ,d$(x,y)  ^  ^du{y)]  j_f  N  f -u‘(x)  ifx<=£>  /0 

J3D  L  “  dn(y)  *U’^dft(2)J  ^  ~  {  uHx)  ifxeR3\£  ('3) 

where  da(y)  is  the  surface  measure  on  dD.  Substituting  (2.2)  in  (2.3)  and  using  (1.8)  we  have: 

1  f  T  d(e~ik{^])  du(y)  ....  ' 

Fo(i,k,a)  =  ±  u(y)"TTl — "  -  Fne",ki'  (2.4) 

4*  JdD  [  -  dn(y)  dn(y) 

Let  B  —  {x  €  R3  |  ||x||  <  1},  dB  be  the  boundary  of  B  and  d\  be  the  surface  measure  on  dB , 
we  denote  with  L'(dB.dX)  the  space  of  square  integrable  complex  functions  with  respect  to  the 
measure  d\.  For  g(x)  €  L2(dB,d\)  and  (2.4)  interchanging  the  integrals  we  have: 


1  f  \  ,  ,dv(y)  du{y)— 
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where 


v(y)  =  [  g{x)eik&'*>d\(x)  (2.6) 

JdB 

It  is  easy  to  see  differentiating  under  the  integral  sign  that  v(y)  is  a  solution  of  the  Helmholtz 
equation  for  every  y  6  R3. 

Let  Ei ,  E-2,  S3  be  the  sets  of  the  eigenvalues  of  the  Laplace  operator  inside  the  domain  D  with  the 
Dirichlet  boundary  condition  (1.5),  the  Neumann  boundary  condition  (1.6)  or  the  mixed  boundary 
condition  (1.7)  respectively.  We  give  the  following  definitions: 

Definition  2. 1  Given  - k 2  £  Ei,  let  uq(y)  be  the  unique  solution  of  the  following  boundary  value 
problem: 

(A  +  k2)tvi(y)  =  0  yeD  (2.7) 


m{y)  =  — 


»*llyll 


y  e  dD 


We  say  that  D  is  an  Herglotz  domain  with  respect  to  the  Dirichlet  boundary  condition  if  there 
exists  <?i(x)  6  L-(dB,d\)  such  that 


w\(y)  =  /  ffi(i)e,'fc(1^)dA(i) 
JdB 


In  a  similar  way  we  define: 

Definition  2.2  Given  -kr  £  Si,  let  um(j/)  be  the  unique  solution  of  the  following  boundary  value 
problem: 

(A  +  *2)u;2(y)  =  0  y€D  (2.10) 

dw^(y)  d  {  e~ik^\ 

- x  v  I  1  a  r\  ( r\  1  1  \ 


y  6  dD 


(2.11) 


dn(y)  dn(y)  ^  *||y|| 

We  say  that  D  is  an  Herglotz  domain  with  respect  to  the  Neumann  boundary  condition  if  there 
exists  y2(x)  6  L 2(dB,d\)  such  that 


«2(y)=  f  92(x)e'k(^d\{x) 
JdB 


(2.12) 


Definition  2.3  Given  \  €  R  and  —  k~  £  S3,  let  ^(y)  be  the  unique  solution  of  the  following 
boundary  value  problem: 


dw^fy)  d 


yeD 

(2.13) 

kM  / 

y  €  dD 

(2.14) 

We  say  that  D  is  an  Herglotz  domain  with  respect  to  the  mixed  boundary  condition  if  there  exists 
g3(i)  €  L2(7  3,  dX)  such  that 


w3(y)  =  [  g3(i)eik^d\(i)  (2.15) 

JdB 

To  our  knowledge  it  is  not  known  a  characterization  of  the  Herglotz  domains,  however  it  is  easy  to 
see  that  the  class  of  Herglotz  domains  is  not  empty.  In  fact  a  straightforward  computation  shows 
that  the  sphere  of  center  the  origin  is  an  Herglotz  domain  in  the  sense  of  Definition  2.1,  2.2,  2.3, 
moreover  the  numerical  experience  of  section  5  can  be  regarded  as  experimental  evidence  that  the 
domains  considered  satisfy  the  previous  definitions. 

Since  now  on  we  consider  only  domains  D  that  satisfy  the  appropriate  Herglotz  condition  that  is 
Definition  2.1  or  2.2  or  2.3.  In  the  case  of  the  inverse  problem  for  the  acoustically  soft  obstacle 
from  (2.5)  using  (1.5)  and  (2.S)  we  have: 


[  F0(x,  =  y 

JdB  * 

Vo6  dB 

(2.16) 

Reasoning  in  the  same  way  we  have: 

/  F0(i,k,a)g2(i)d\(i)  =  y 

JdB  K 

V  a  G  dB 

(2.17) 

for  the  acoustically  hard  obstacle  and 

[  F0ii.k.a)g3(i)d\(i)  =  y 
JdB  * 

dB 

'2.18) 

for  the  obstacle  characterized  by  an  acoustic  impedance. 

The  numerical  method  for  the  inverse  problem  for  the  acoustically  soft  obstacle  is  based  on  the 
relations  (2.16),  (2.9),  (2.8)  that  connect  the  data  that  is  the  far  fields  to  the  unknown  dD.  In  a 
similar  way  we  will  exploit  (2.17),  (2.12),  (2.11)  to  solve  the  inverse  problem  for  the  acoustically 


hard  obstacle  and  (2.18),  (2.15),  (2.14)  to  solve  the  inverse  problem  for  the  obstacle  characterized 
by  an  acoustic  impedance. 

3-  The  numerical  method 

Given  D  C  R3  and  the  boundary  condition  (1.5)  or\  1.6)  or  ( 1.7)  satisfied  by  u  on  dD  we  will  exploit 
numerically  the  analytic  relations  derived  in  section  2  as  follows.  Since  most  of  the  content  of  this 
section  is  independent  of  the  boundary  conditions  chosen,  in  order  to  fix  the  ideas,  we  consider 


the  inverse  problem  for  acoustically  soft  obstacles  that  is  the  relations  (2.16),  (2.9),  (2.3),  when 
necessary  we  will  comment  on  the  peculiar  features  of  the  corresponding  problems  for  acoustically 
hard  obstacles  or  obstacles  characterized  by  an  acoustic  impedance.  Cur  general  strategy  car.  be 
summarized  in  three  points: 

(i)  use  (2.16)  to  go  from  the  knowledge  of  the  fa.  ~  Ids  Fq  to  the  Hergiotz  kernel  gi(i)  of  the 
domain  D 

(ii)  use  (2.9)  to  go  from  the  knowledge  of  the  Hergiotz  kernel  gi(i)  to  the  Hergiotz  wave  function 
wi(y) 

(iii)  use  (2.8)  to  go  from  the  knowledge  of  the  Hergiotz  wave  function  wl(y)  to  the  boundary  of 
the  obstacle  OD 

More  precisely  given  D,  let  =  (o,  €  dB  |  i  ~  1, 2,  •  •  • ,  iV}  be  the  set  of  directions  of  the  incoming 
waves,  fl3  =  {fcj  6  R  |  -  k;  £  Et  ,  i  =  1,2,  be  the  set  of  the  non-resonant  wave  numbers 

of  the  incoming  waves  and  fl3  =  {i,  6  dB  |  i  =  1, 2,  •  •  • ,  M)  be  the  set  of  directions  where  the  far 
fields  Fo  are  measured.  For  i  =  1,2,3  we  assume  that  the  elements  of  fl,  are  distinct.  The  data  of 
our  inverse  problems  will  be  the  numbers  P,.;,  j,  i  =  1,2,  •  •  •,  N,  t'i  =  1,2,  •  •  •,  jV1(  j  =  1,2,  •  •  • ,  M 
that  represent  the  measurements  of  Fo(i;, £,,, a,).  In  the  numerical  experience  shown  in  section  5 
these  data  are  obtained  solving  numerically  the  direct  problem  (1.3),  (1-4),  (1.5). 

Let  (6,  <p)  be  the  polar  angles  so  that 

i  =  i{6,  o)  =  (sin#cosd>,  sin#sind>,  cosd)T  (3-1) 

and  Ulm(i)  =  ~/im  P™ (cos  0)  cos  mo.  Vim(x)  -  -yimP[n(cos0)  s'mmp,  /  =  0, 1,  2,  -  -  -  ,  m  —  0,1,---,/ 
be  the  spherical  harmonics,  where  are  the  normalization  factors  in  L2(dB,d\)  and  P/7*  are  the 
Legendre  functions.  It  is  well  known  that  {Uim,Vim}  1=0.1.  •  is  an  orthonormal  complete  set  of 

m  =0, 1 .  •  •  ,J 

V(dB,dX). 

From  the  data  our  computation  proceeds  as  follows: 

Step  1.  From  the  data  to  the  coefficients  of  the  expansion  in  spherical  harmonics  of  the  far  field 
F0(x, *<,,«,)  ,  ii  =  1,2,-  -  - , - V !  ,  t  =  1,2,  "-..V. 

Let  us  first  consider  the  case  fl3  =  {z,  |  i  =  1,2,  •  •  •,  M}.  Given  an  integer  Lmai  >  Owe  assume  that 
the  far  field  P0(i,  kil,ai)  can  be  approximated  by  a  truncated  expansion  in  spherical  harmonics. 


that  is: 


Lmat  1  I 

E  E  fSL.i  Utm(i)  +  Y,  E  F0j:m.2V>m(i) 

1=0  m=0  i=l  m=  1 
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we  note  that  the  unknown  coefficients  appearing  in  the  truncated  expansion  (3.2)  are  ( Lmax  +  l)2 
complex  constants.  To  determine  these  coefficients  we  impose  that 


— ' rn a  s  •  “nio* 

£  £c £  £  f&LjlW*,)  -  Km  .  i-i.V 


/=0  m=0 


i=l  m=l 


The  equations  (3.3)  are  a  linear  system  of  M  equations  in  ( Lmax  4-  l)2  unknowns,  so  that  in  order 
to  determine  the  unknowns  we  need  M  >  ( Lmax  +  l)2.  The  condition  number  of  the  linear  system 
(3.3)  grows  dramatically  with  the  number  of  equations  M,  so  that  its  numerical  solution  become 
practically  unfeasible  even  for  small  values  of  Zmax  such  as  Lma x  —  4,6,8. 

The  use  of  regularization  procedures  such  as  Tichonov  regularization  to  solve  (3.3)  are  helpful  but 
insufficient  to  cure  the  ill  conditioning.  We  note  that  the  right  hand  side  of  (3.3)  is  supposed  to 
represent  actual  physical  measurements  that  are  affected  by  significant  experimental  errors  so  that 
ill  conditioning  is  a  true  challenge.  This  ill  conditioning  problem  is  solved  reducing  the  number  of 
unknowns  to  be  determined  and  reducing  in  a  similar  way  the  number  of  equations  to  be  solved. 
This  is  obtained  exploiting  the  special  features  of  the  linear  system  (3.3). 


Step  2.  From  the  far  fields  Fo{i, ,a;)  i\  =  1, 2,  •••, Ni  i  =  1, 2,  ■  •  • , N  to  the  Herglotz  kernel 

Let  <?,,(x)  €  L2(dB,d\)  be  the  Herglotz  kernel  associated  to  the  boundary  condition  (1.5),  the 
domain  D  and  the  wave  number  ku  €  Qi-  We  assume  that  gi,(i)  can  be  appoximated  by  a 
truncated  expansion  in  spherical  harmonics,  that  is: 

L7  i  L1  i 

=  L  E  9im +  E  E  (3.4) 

(=0  m=0  /=  1  1 

where  0  <  L3  <  Lmax  is  an  integer.  Using  the  orthogonality  pr  .<-’s  of  the  spherical  harmonics 
the  relation  (2.16)  will  be  approximated  with: 


£  £ 


1=0  m=G 


n 1 1  —  _ 

y.yim2  l 
K'\ 


i  =  1,2  V 


The  equations  (3.5)  are  a  linear  system  of  N  equations  in  the  ( Lg  +  l)2  unknown  coefficients 
{lyfmi’l/imz}’  so  t^at  *n  order  to  determine  the  unknowns  we  need  N  >  ( Lg  +  l)2.  We  remark  that 
the  unknown  coefficients  {<7fmp  are  complex  numbers.  The  number  of  incoming  waves  N  can 
be  drastically  reduced  if  the  unknown  coefficients  to  be  determined  are  reduced  in  a  corresponding 
manner.  This  reduction  can  be  achieved  if  we  make  some  a  priori  assumptions  about  the  symmetries 


of  the  obstacle  D ,  and  we  assume  than  the  same  symmetry  is  conserved  in  the  Herglotz  kernel 
In  this  paper  we  consider  only  two  choices: 

(i)  dD  is  cylindrically  symmetric  with  respect  the  z-axis,  that  is  dD  =  {x  =  f(0)x{9,d>)  |  0  < 
9  <  z  ,  0  <  <p  <  2x}  for  some  function  /.  In  this  case  if  we  assume  the  same  symmetry  for 
gix{i)  we  have  =  0  if  m  >  0  and  g }^2  =  0  for  every  l,m  so  that  the  linear  system  (3.5) 
has  only  Lg  +  1  unknowns  and  as  a  consequence  we  need  only  N  >  Lg  4-  1  incoming  waves 

(ii)  dD  is  cilindrically  symmetric  with  respect  to  the  z-axis  and  symmetric  with  respect  to  the 


equator.  Arguing  as  in  (i)  we  can  conclude  that  the  linear  system  (3.5)  has  only 


+  1 


unknowns.  With 


we  have  denoted  the  integer  part  of  . 

When  it  is  necessary  the  system  (3.5)  is  solved  using  a  regularization  procedure.  We  note  that 
when  it  is  assumed  some  symmetry  about  dD  and  these  symmetries  are  exploited  as  previously 
shown  a  large  number  of  the  coefficients  {/oY.m.p  -^o.Vm  z)  °f  t^le  system  (3-5)  determined  in  Step 
1  are  multiplied  by  unknowns  that  are  assumed  zero. 

Step  3.  From  the  Herglotz  kernel  gix(i)  to  the  Herglotz  wave  function  uqti,(y),  i\  =  1,2 
From  (2.9)  and  (3.4)  an  explicit  computation  gives  us  an  approximated  expression  of  U7it;,(y),  that 


is: 


wi,ii(y)  =  4?r  l 


Lg  l  Lg  l 

E  L  iM  w«(i>  +  E  E  iis/iD'^® ;  (3.6) 


/=0  m=0  1=  1  m=  1 

where  y  =  y/||y||  and  ji  is  the  spherical  Bessel  function  of  order  /. 


Step  4.  From  the  Herglotz  functions  uiijM(y),  i'i  =  1,2,  to  the  boundary  of  the  obstacle 

dD. 

For  simplicity  we  assume  that  there  exists  0  <  a  <  b  <  +oo  and  a  smooth  function  such 

that  a  <  <  b,  V9,<t>  and  we  have  dD  =  {( r,9,<p )  \  r  =  0  <  9  <tt,  0  <  <t>  <  2t}. 

The  reconstruction  of  the  boundary  dD  from  the  Herglotz  functions  is  based  on  the  relation  (2.3)  in 
the  case  of  the  acoustically  soft  obstacle,  or  the  relation  (2.11)  in  the  case  of  the  acoustically  hard 
obstacle,  or  the  relation  (2.14)  in  the  case  of  the  obstacle  characterized  by  an  acoustic  impedance. 
The  relation  (2.8)  must  be  interpreted  as  a  nonlinear  equation  that  defines  implicitely  /  as  a 
function  of  9  and  4>,  the  relations  (2.11),  (2.14),  due  to  the  presence  of  the  d/dn  term,  are  nonlinear 
expressions  involving  /,  |^,|^  that  is  first  order  partial  differential  equations.  So  that  we  expect 
that  the  reconstruction  of  dD  in  the  acoustically  soft  case  should  be  easier  than  in  the  remaining 
cases.  In  section  4  we  will  show  an  ad  hoc  procedure  to  exploit  (2.8),  here  we  restrict  our  attention 


to  a  general  procedure  to  obtain  dD  that  can  be  applied  always.  We  approximate  dD  with  a 
truncated  series  of  spherical  harmonics: 

l  Lp  l 

0  =  (3.7) 

(=0  m=0  1=1  m=l 

where  Lp  >  0  is  chosen  depending  on  dD ,  i.e.  simple  obstacles  can  be  reconstructed  with  small 
I?,  that  is  Lp  =  4,6.  Moreover  if  it  is  assumed  that  dD  has  some  of  the  symmetries  previously 
considered  these  symmetries  can  be  easily  translated  into  properties  of  the  coefficients  {c(mi,C(m2}, 
that  is  the  appropriate  coefficients  can  be  chosen  zero.  Let  c  =  {cfmi,c/m2},  0  <  /  <  Lp,  0  <  m  <  l 
be  the  ( Lp  +  l)2  dimensional  real  vector  of  the  unknown  coefficients,  the  unknown  boundary  is 
obtained  minimizing  with  respect  to  c  the  following  functions 

A  1  f2*  r  e~'k'J  2 

A(c)  =  V  t-  /  d<t>  sin0  wUl  +  -7—7-  dO  (3-8) 

i1=i  Kii  Jo  KfiJ 

when  the  acoustically  soft  obstacle  is  considered, 

A  1  f2*  r  d  e~ik'\!  2 

^)=Ej-yo  Tn(W2-i'+—)  di  (3-9) 

when  the  acoustically  hard  obstacle  is  considered, 

A  1  f2n  r  d  e~ikiif  2 

h (£)=Y2~J0  d<i> J0  s[nd  (l  +  Xg^)(w3,il  + dQ  (3.10) 

when  the  obstacle  with  acoustic  impedance  is  considered. 

In  (3.8),  (3.9),  (3.10)  the  Herglotz  wave  functions  obtained  in  Step  3  are  computed  in  (r  = 
/(0,0),0, 0),  ^  is  the  normal  derivative  with  respect  to  the  surface  r  =  )  and  /  is  ap¬ 

proximated  with  (3.7).  The  factor  in  f3.8),  (3.9),  (3.10)  is  chosen  to  make  the  addenda  corre- 
sponding  to  different  values  of  k  of  the  same  order  of  magnitude,  finally  the  integrals  appearing  in 
(3.8),  (3.9),  (3.10)  are  approximated  by  some  simple  quadrature  rule.  We  observe  that  the  functions 
(c) ,  i  =  1,2,3  are  non  negative  functions  and  that,  if  we  neglect  the  effects  of  the  approximations 
introduced,  the  surface  r  =  /(#.©)  corresponds  to  a  point  where  f,(c)  i  =  1,2,3  is  zero,  that  is 
r  =  is  a  global  minimizer  of  /,(c).  We  note  that  in  general  the  surface  dD  is  only  a  proper 

subset  of  the  set  of  points  satisfying  the  relations  (2.8)  or  (2.11)  or  (2.14). 

When  the  minimization  of  /,(c)  does  not  give  a  satisfactory  reconstruction  of  dD  we  minimize  Pi(c) 
instead  of  /j(c),  where: 

P(c)  =  rt(c)  +  q(c)  i  =  1.2,3  (3.11) 
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and  q(c )  is  a  penalization  term.  The  penalization  term  q(c )  is  chosen  as  follows 


Til  Tin 


=  PiAWi’Vj)-#)2 

:  =  1  j=  1 


(3.12) 


where  >  0  are  weight  factors,  /*■  are  approximated  values  of  /  in  the  direction  (0, ,  6,)  that 
we  assume  known.  For  the  acoustically  soft  obstacles  the  values  f*-  can  be  obtained  solving  (2.8) 
with  9  —  9{,  <p  =  <pj  as  a  nonlinear  equation  for  /.  In  the  remaining  cases  the  nonlinear  partial 
differential  equations  (2.11),  (2.14)  become  nonlinear  equations  when 


df 

89 


(3.13) 


For  a  large  class  of  surfaces  (3.13)  is  satisfied  in  some  special  points  such  as  the  North  Pole  (i.e.  9  — 
0,  4>  arbitrary),  the  South  Pole  (i.e.  9  =  t,  <p  arbitrary)  or  the  Equator  (i.e.  9  =  ~ ,  <p  arbitrary). 
At  these  points  (2.11),  (2.14)  can  be  solved  as  nonlinear  equations  to  obtain  approximate  values 
/*  of  /  to  be  used  in  the  penalization  term  (3.12). 

The  minimization  of  /{(c)  or  P,(c)  is  performed  with  one  of  the  following  three  algorithms:  DB- 
CONF  [17],  that  is  a  quasi  Newton  local  minimization  algorithm,  SIGMA  [18], [19],  that  is  a  global 
minimization  stochastic  algorithm,  or  DUNLSJ  [20],  that  is  a  non  linear  least  squares  algorithm. 
The  quasi-Newton  algorithm  performs  only  a  local  minimization  but  is  computationally  cheaper 
than  the  global  minimization  alghorithm.  The  nonlinear  least  squares  algorithm  is  specifically 
suited  for  the  minimization  of  /.(c)  i  =  1,2.3. 


4-  The  reconstruction  of  OD  in  the  acoustically  soft  case 

In  section  3  we  have  observed  that  (2.8)  can  be  interpreted  as  a  nonlinear  equation  defining  / 
implicitely  as  a  function  of  9  and  o.  So  that  differentiating  (2.8)  with  respect  to  9  and  0  we  have: 


dw\  ,.  ...  1  .  e~,ku  ' 


0j_  duti.,,  _ 

09  89 


=  0  ,  0  <  9  <  ft,  O<0<2tr 


8£  flwt.,,  _ 

8<p  do 


,  O<0<~,  0  <  <p  <  2t 


(4.1) 


(4.2) 


8r  '  '  *„  '  P 

Equations  (4.1),  (4.2)  can  be  interpreted  as  ordinary  differential  equations  for  the  unknown  /. 
We  note  that,  in  the  reconstruction  procedure  considered  here,  /  is  not  approximated  with  a 
truncated  series  in  spherical  harmonics  as  in  (3.7)  and  that  the  choice  of  several  values  of  k  (i.e. 
fcM,  i!  =  1,2,  --,.V1)  does  not  play  any  role.  In  order  to  obtain  f(6,<j>)  we  proceed  as  follows: 
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(i)  for  fixed  9  and  d>,  let  say  9  =  <t>  =  0,  we  solve  (2.8)  as  a  nonlinear  equation.  Let  /0 0  be  the 
solution  found 

(ii)  given  /oo  we  solve,  for  0  <  9  <  x,  the  differential  equation  obtained  by  taking  the  real  part  of 
(4.1)  with  the  initial  condition 

/(0,0)  =  /oo  (4.3) 

Let  f0(9, 0)  be  the  solution  found.  We  verify  that  fo(9,Q)  satisfies  (2.8) 

(iii)  given  /o(0,O)  we  solve,  for  0  <  6  <  2x,  the  differential  equation  obtained  by  taking  the  real 
part  of  (4.2)  with  the  initial  condition 

/(M)  =  /o(M)  (4.4) 

Finally  we  verify  that  the  obtained  in  this  way  satisfies  (2.8). 

The  differential  problems  considered  in  (ii),  (iii)  are  initial  value  problems  for  scalar  differential 
equations  and  are  solved  numerically  using  a  Runge-Kutta-Fehlberg  method  (i.e.  the  subroutine 
RKF45  [21]). 

The  problem  considered  in  (iii)  is  performed  only  a  finite  number  of  times  corresponding  to  a 
discretization  of  the  variable  9. 

We  note  that  if  the  obstacle  is  cylindrically  symmetric  with  respect  to  the  z  axis  (iii)  is  not  necessary 
and  the  problem  is  solved  after  performing  (i),  (ii).  Moreover  the  differential  equations  of  (iii)  for 
different  values  of  9  are  independent  one  from  the  other  and  can  be  solved  in  parallel.  The  set 
of  solutions  of  (2.8)  in  general  contains  the  surface  f(9,<f>)  as  a  subset  so  that  when  performing 
(i),  (ii),  (iii)  only  trajectories  that  appear  to  define  a  closed  surface  should  be  considered.  Finally 
the  solution  of  the  inverse  problem  based  on  the  differential  equations  (4.1),  (4.2)  is  limited  to  the 
acoustically  soft  obstacles  and  appears  to  be  more  sensible  to  error  in  the  data  than  the  minimization 
procedure  described  in  section  3.  However  its  computational  cost  is  very  small  compared  to  the 
minimization  procedures  previously  described. 

5.  The  numerical  experience 

In  this  section  we  describe  the  numerical  results  obtained  using  the  numerical  methods  described 
in  sections  3  and  4  on  several  test  problems.  When  we  consider  the  boundary  condition  (1.7)  we 
choose  x  =  1.  The  surfaces  dD  considered  are  the  follo%ving  ones: 

(^x)'  +  (|i/)2  +  =  1  Oblate  Ellipsoid  (5.1) 
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where 


((fx)!  +  (fj)3)s  +  x'°  =  i 

Short  Cylinder 

(5.3) 

(X2  +  i-f  +  (fx)'°  =  1 

Long  Cylinder 

(5.4) 

2  .  ,  .  1  .  i  ..  i 

r  =  -(cos*  9  +  -  sin"  9)? 

3  4  ’ 

Vogel's  Peanut 

(5.5) 

r  =  1  —  i  cos  2 9 

2 

Horizontal  Platelet 

(5.6) 

r  =  1  -j-  ^  cos  29 

2 

Vertical  Peanut 

(5.7) 

3. 17  i 

r  =  -(—  +  2  cos  30)3 

Pseudo  Apollo 

(5.8) 

Corrugated  Ellipsoid 

(5.9) 

=(0!+ (3«*>,r* 

Ellipsoid 

(5.10) 

/,2sm0,io 

V3  ) 

Corrugated  Cylinder 

(5.11) 

r  =  ^  +  H(4>)  sin2  6 

Lt 

Corrugated  Platelet 

(5.12) 

h(<j>)  =  ((^  cos<£)2  4-  sin2  <t>)  % 

=  (R^  +  Ah  cos  4<p  +  Bh  cos  8<t>  +  Ch  cos  16<£)2 


(5.13) 

(5.14) 


and 

Ah  =  rkfh  ;  =  :  Ch  =  ~TT4fh  ]  ^  =  i-(^  +  5*+cv)  (5-i5) 

fh  is  the  corrugation  parameter.  In  our  numerical  experience  fh  =  0.2.  The  obstacles  D  corre¬ 
sponding  to  (5.1),  (5.2),  (5.3),  (5.4)  are  convex  bodies  symmetric  with  respect  to  the  z-axis  and 
to  the  equator,  the  obstacles  corresponding  to  (5.5),  (5.6),  (5.7),  (5.8)  are  non  convex  but  they 
mantain  the  symmetry  with  respect  to  the  z-axis.  Finally  the  obstacles  D  corresponding  to  (5.9), 
(5.10),  (5.11),  (5.12)  in  general  are  non  convex  and  non  symmetric  with  respect  to  the  z-axis.  We 
observe  that  a  characteristic  length  L  of  the  obstacles  can  be  chosen  to  be  one.  The  Tables  5.1,  5.2, 
5.3,  5.4  show  some  numerical  results  obtained  with  the  methods  described  in  the  previous  sections. 
In  the  Eh,  and  Emax  columns  of  those  tables  we  use  the  notation  a(/?)  to  mean  a  ■  10^.  Table  5.1 
summarizes  the  results  obtained  with  the  obstacles  (5.1),  (5.2),  (5.3),  (5.4),  Table  5.2  summarizes 
the  results  obtained  with  the  obstacles  (5.5),  (5.6),  (5.7),  (5.8),  finally  Table  5.3  summarizes  the 
results  obtained  with  the  obstacles  (5.9),  (5.10),  (5.11),  (5.12). 
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We  remark  that  in  the  reconstructions  presented  in  Table  5.3  the  obstacles  are  considered  as  general 
surfaces,  that  is  there  is  no  use  of  symmetries  in  the  reconstruction  procedure.  In  Figs.  5.1,  5.2... .,5. 9 
the  reconstructions  denoted  with  (*)  in  the  Tables  5.1,  5.2,  5.3  are  shown. 

Let  —  (ijjg,  i  =  1,2...., 35  ,  j  =  0, 1 , . . . , 3 5  and  let  /(#,,  q:)  be  the  exact  values  of  the 

surfaces  given  by  (5.1),  (5.2), ...,(5.12)  and  <pj)  be  the  corresponding  values  obtained  by  the 
reconstruction  procedure  of  sections  3  and  4.  In  the  Tables  5.1,  5.2,  5.3,  5.4  we  use  as  performance 
index  the  relative  L2  error,  that  is: 


Er,  = 


(/(0,0)  -  /c(0,0))2  +  (/(*,  0)  -  /e(x,  0))2  +  £?-i  Ef=0 -  m,*:))2 


1/2 


/( 0,0)2  +  /(t,0)2  +  M'tj) 


1/2 


(5.16) 


Table  5.1  Axially  symmetric  convex  obstacles 
f  max  —  Lg  —  8  fix  —  A\  J 3  “  /It 


Object 

Recon¬ 

struction 

Boundary 

condition 

k\ 

^2 

Reconstruction 

method 

Penalizati 
North  Pole 

on  term 
Equator 

e 

Ei,- 

Oblate  Ellipsoid 

1 

Dirichlet 

- 

3 

- 

RKF45 

- 

- 

0.0 

0.868(-4) 

n 

2 

77 

- 

3 

- 

RKF45 

- 

- 

0.05 

0.703(-3) 

77 

3(*) 

Neumann 

4 

2 

no 

DBCONF 

no 

no 

0.0 

0.546(-2) 

77 

4(*) 

77 

4 

2 

no 

DBCONF 

no 

no 

0.05 

0.880(-2) 

77 

5 

Mixed 

4 

2 

no 

DBCONF 

no 

no 

0.0 

0.365(-2) 

77 

6 

77 

4 

2 

no 

DBCONF 

no 

no 

0.05 

0.589(-2) 

Prolate  Ellipsoid 

7 

Dirichlet 

- 

3 

- 

RKF45 

- 

- 

0.0 

0.256(-4) 

■7 

8 

77 

- 

3 

- 

RKF45 

- 

- 

0.05 

0.966(-2) 

77 

9 

Neumann 

4 

2 

no 

DBCONF 

no 

no 

0.0 

0.172(-1) 

77 

10 

77 

4 

2 

no 

DBCONF 

no 

no 

0.05 

0.157(-1) 

,7 

11 

Mixed 

4 

2 

no 

DBCONF 

no 

no 

0.0 

0.432(-2) 

77 

12 

*7 

4 

2 

no 

DBCONF 

no 

no 

0.05 

0.869(-2) 

Short  Cylinder 

13 

Dirichlet 

- 

3 

- 

RKF45 

- 

- 

0.0 

0.629(-2) 

77 

14 

77 

- 

3 

- 

RKF45 

- 

- 

0.05 

0.147(-1) 

77 

15 

Neumann 

6 

2 

no 

SIGMA 

no 

no 

0.0 

0.304(-l) 

♦7 

16 

•7 

6 

2 

no 

SIGMA 

yes 

yes 

0.05 

0.626(-l) 

77 

17(*) 

Mixed 

6 

2 

no 

SIGMA 

yes 

yes 

0.0 

0.306(-l) 

77 

isn 

77 

6 

2 

nq 

SIGMA 

yes 

yes 

0.05 

0.468(-l) 

Long  Cylinder 

19(+) 

Dirichlet 

- 

3 

- 

RKF45 

- 

- 

0.0 

0.302(-2) 

77 

20(*) 

77 

— 

3 

- 

RKF45 

- 

- 

0.05 

0.369(-2] 

77 

21 

Neumann 

6 

2 

3 

SIGMA 

yes 

yes 

0.0 

0.543(-l) 

77 

22 

77 

6 

2 

3 

SIGMA 

yes 

yes 

0.05 

0.482(-l) 

77 

23 

Mixed 

6 

2 

no 

SIGMA 

yes 

yes 

0.0 

0.415C-1) 

77 

24 

77 

6 

2 

no 

SIGMA 

yes 

yes 

0.05 

0.425(  1) 

or  the  relative  L°°  error,  that  is: 


_  max{!/(0, 0)  -  /e(0,0)l,j/(x,0)  -  /e(x,  Q)|,  \f(9it  o:)  -  } 


E-max  — 


Let  us  define  the  sets: 


max{ |/(0, 0)|, |/(x, 0)|, !/(#,, <P>)|  ,  > 


Ai  =  {(<?;, 0)  |  9i  =  ,  i  =  0,1,---,  10} 


.  X 


a2  =  {(0.-, 4>j)  i  =  =  i  i  =  o, l, •••,8} u  {o,o} u  {x,o} 

A3  =  {(<?!,<£;,)  |  Ox  -  ,i=  7  ;  <pj  =  ,  j  =  0, 1, •  •  • , 6}  U  {0,0}  U  {x,0} 


(5.17) 

(5.13) 

(5.19) 

(5.20) 


Table  5.2  Axially  symmetric  non-convex  obstacles 
Imar  “  ^ g  ~  8  ,  O 1  —  A \  ,  D 3  ”  Ao 


Object 

Recon¬ 

struction 

Boundary 

condition 

1 

*2 

Penalizati 
North  Pole 

on  term 
Equator 

e 

eL3 

Vogel’s  Peanut 

in 

Dirichlet 

4 

3 

nd 

SIGMA 

no 

no 

0.0 

0.128(-1) 

11 

2(*) 

rt 

4 

3 

nc 

SIGMA 

no 

no 

0.05 

0.175(A) 

n 

3 

Neumann 

4 

2 

no 

SIGMA 

no 

no 

0.0 

0.596(A) 

ii 

4 

’1 

4 

2 

no 

SIGMA 

no 

no 

0.05 

0.176(A) 

ii 

5 

Mixed 

4 

2 

nc 

SIGMA 

no 

no 

0.0 

0.273(A) 

ii 

6 

•1 

4 

2 

ss 

SIGMA 

no 

no 

0.05 

0.128(A) 

Horizontal  Platelet 

1 

Dirichlet 

4 

nr 

no 

SIGMA 

no 

no 

11 

8 

ii 

4 

3 

nc 

SIGMA 

no 

no 

failure 

11 

9 

Neumann 

4 

2 

nc 

SIGMA 

no 

no 

1 1 

10 

11 

4 

2 

nc 

SIGMA 

no 

no 

11 

nn 

Mixed 

4 

2 

nc 

SIGMA 

no 

no 

m 

11 

12  n 

4 

2 

no 

SIGMA 

no 

no 

Vertical  Peanut 

tmsm 

Dirichlet 

4 

3 

SIGMA 

yes 

yes 

0.0 

0.172(A) 

11 

mm 

4 

3 

TO 

SIGMA 

yes 

yes 

0.05 

failure 

11 

15 

Neumann 

4 

2 

3 

SIGMA 

no 

no 

0.0 

0.725(A) 

11 

16 

4 

2 

3 

SIGMA 

no 

no 

0.05 

failure 

11 

17 

Mixed 

4 

2 

3 

SIGMA 

no 

no 

0.0 

0.649(A) 

11 

18 

'1 

4 

2 

3 

SIGMA 

no 

no 

0.05 

failure 

Pseudo  Apollo 

19 

Dirichlet 

4 

3 

nc 

SIGMA 

no 

no 

0.0 

0.602(A) 

ii 

20 

4 

3 

no 

SIGMA 

no 

no 

0.05 

0.114(4-0) 

11 

21(*) 

Neumann 

4 

2 

3 

SIGMA 

yes 

no 

0.0 

0.359(A) 

11 

22(*) 

T1 

4 

2 

3 

SIGMA 

yes 

no 

0.05 

0.325(A) 

11 

23 

Mixed 

4 

2 

3 

SIGMA 

yes 

no 

0.0 

0.551(A) 

11 

24 

4 

2 

3 

SIGMA 

yes 

no 

0.05 

0.32  If  - 1 ) 
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In  our  numerical  experience  we  choose 

Qi  =  At  (5.21) 

when  the  obstacles  of  Tables  5.1,  5.2  are  considered  or 

ftt  =  A3  (5.22) 

when  the  obstacles  of  Table  5.3  are  considered.  T1  et  Qo  is  either  {2},  {3},  or  {2,3}  (see  Tables). 
Finally  when  the  obstacles  of  Tables  5.1,  5.2  considered  fl3  is  given  by  Ao  and  when  the  obstacles 
of  Tables  5.3  are  considered  03  is  given  by  A3.  We  observe  that  with  these  choices  the  resonance 
condition  (1.9)  is  satisfied.  The  far  field  data  corresponding  to  these  choices  are  obtained  by  solving 
numerically  the  corresponding  direct  problems  that  is  the  boundary  value  problems  (1.3),  (1.4), 
(1.5)  or  (1.3),  (1.4),  (1.6)  or  (1.3),  (1.4),  (1.7)  using  a  T-matrix  approach  [16]. 


Table  5.3  Generic  obstacles 


Reconstruction  method  =  ”D UNLSJ”,  without  penalization  term, 
Lmax  “  -^7  ~  6}  Lp  =  4;  k\  —  2,  Qi  —  A3,  fl3  =  A3 


Object 

ReconstructionjBoundary  condition 

e 

Ev 

Corrugated  Ellipsoid 

1 

Dirichlet 

0.0 

0.276(-l) 

11 

2 

11 

0.02 

0.428(-l) 

11 

3(*) 

Neumann 

0.0 

0.287(-l) 

” 

4(*) 

71 

0.02 

0.417(-1) 

17 

5 

Mixed 

0.0 

0.325(-l) 

11 

6 

11 

0.02 

0.133(4-0) 

Ellipsoid 

f 

Dirichlet 

0.0 

0.291(-1) 

M 

3 

'1 

0.02 

0.296(-l) 

?7 

9 

Neumann 

0.0 

0.303(-l) 

11 

10 

0.02 

0.308(-l) 

u 

11 

Mixed 

0.0 

0.152(+0) 

n 

12 

'1 

0.02 

0.245(-l) 

Corrugated  Cylinder 

13 

Dirichlet 

0.0 

0.534(-l) 

17 

14 

11 

0.02 

0.652(-l) 

17 

15 

Neumann 

0.0 

0.502(-l) 

'7 

16 

' 1 

0.02 

0.11S(+0) 

11 

!"(*) 

Mixed 

0.0 

0.569(-l) 

11 

13(  *) 

11 

0.02 

0.675(-l) 

Corrugated  Platelet 

ion 

Dirichlet 

0.0 

0.663(1) 

20(  * ) 

11 

0.02 

0.129(4-0) 

? 1 

21 

Neumann 

0.0 

0. 105(4-0) 

11 

22 

0.02 

failure 

'1 

23 

Mixed 

0.0 

0.65  1(-1) 

n 

24 

<1 

0.02 

0.175(4-0) 

33 


Table  5.4  Performance  as  a  function  of  e 
Object:  Long  Cylinder;  Boundary  condition:  Dirichlet; 

Lmax  =  Lg  =  8;  Lp  =  6;  k\  =3;  Qi  =  At  fl3  =  A3 
Reconstruction  method=  ’’SIGMA”  without  penalization  term. 


Reconstruction!  e 

L’mar 

El* 

1 

0.0 

0.611(-1) 

0.286(-l) 

2 

0.05 

0.614(-1) 

0.294(-l) 

3 

0.10 

0.587(-l) 

0.298(-l) 

4 

0.20  0.525(-l) 

0.398(-l) 

5 

0.30  0.865(-l) 

0.548(-l) 

6  j0.40j  0.109(4-0) 

0.655(-l) 

7  0.50  0.123(4-0) 

0.714(-1) 

8  0.60  0.131(4-0) 

0.752(-l) 

To  the  data  =  Fo(i;,  is  added  a  random  error  term,  that  is  Fi,iuj  is  substituted  with 


F i.ijj  —  Ft,tuj  4- 


(5.23) 


where  e  >  0  is  a  parameter  and  (,'  is  a  random  number  uniformly  distributed  in  [—1, 1]. 

In  our  numerical  experience  we  have  Lmax  =  L3  =  6  or  8,  Lp  =  4  or  6  (see  Tables).  Finally  in 
Table  5.4  we  show  the  performance  of  our  algorithms  for  increasing  values  of  e  in  the  case  of  the 
acoustically  soft  long  cylinder.  The  method  based  on  the  global  minimization  algorithm  SIGMA 
appears  to  be  the  most  powerful  one  at  the  price  of  higher  computational  cost.  The  computations 
previously  described  have  been  performed  on  a  VAX  6310  with  VMS  Operating  System. 
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ORIGINAL 


RECONSTRUCTION  19 
of  table  S.l 


RECONSTRUCTION  20 
of  table  5 . 1 


Fig.  5.1  Long  Cylinder 


ORIGINAL 


RECONSTRUCTION  3 
of  table  5.1 


RECONSTRUCTION  4 
of  table  5.1 


Fig.  5.2  Obla.te  Ellipsoid 


ORIGINAL 


RECONSTRUCTION  17 
of  table  5.1 


RECONSTRUCTION  18 
of  table  5.1 


Fig.  5.3  Short  Cylinder 
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ORIGINAL  RECONSTRUCTION  1 

of  table  5.2 


Fig.  5.4  Vogel’s  Peanut 


ORIGINAL 


RECONSTRUCTION  21 
of  table  5.2 


Fig.  5.5  Pseudo  Apollo 


ORIGINAL  RECONSTRUCTION  11 

of  table  5.2 


Fig.  5.6  Horizontal  Platelet 


RECONSTRUCTION  2 
of  table  5.2 


RECONSTRUCTION  22 
of  table  5.2 


RECONSTRUCTION  12 
of  table  5.2 
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ORIGINAL 


RECONSTRUCTION  19 
of  table  5.3 


RECONSTRUCTION  20 
of  tabla  5.3 


Fig.  5.7  Corrugated  Platelet 


ORIGINAL 


RECONSTRUCTION  3 
of  table  5.3 


RECONSTRUCTION  4 
of  table  5.3 


Fig.  5.8  Corrugated  Ellipsoid 
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Fig.  5.9  Corrugated  Cylinder 
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Abstract 

A  numerical  method  for  the  three  dimensional  inverse  acoustic  and  electromagnetic 
time  harmonic  scattering  problem  is  presented.  The  far  field  patterns  of  the  Helmholtz 
or  vector  Helmholtz  equations  generated  by  a  known  wave  incident  on  an  obstacle  D  are 
measured.  These  measurements  are  repeated  for  several  incoming  waves.  From  these 
measurements  the  boundary  dD  of  the  obstacle  is  reconstructed.  The  reconstruction 
procedure  proposed  here  generalizes  the  "Herglotz  function  method”  introduced  by 
Colton  and  Monk  [1]  in  the  acoustic  problem  and  is  effective  in  the  so  called  resonance 
region. 


1.  INTRODUCTION 

Let  R3  be  the  three  dimensional  euclidean  space,  z  =  (r,y,z)  g  RJ  be  a  generic 
vector,  (.,.)  will  denote  the  euclidean  scalar  product  and  j|.||  the  euclidean  norm.  Let 
D  C  R3  be  a  bounded  simply  connected  domain  with  smooth  boundary  dD  that 
contains  the  origin.  Let  u'fz)  be  an  incoming  acoustic  plane  wave,  that  is: 

u'(i)  =  (1.1) 

where  k  >  0  is  the  wave  number  and  a  6  R5  is  a  fixed  unit  vector.  Let  us  denote  with 
u’(x)  the  acoustic  field  scattered  by  the  obstacle  D  and  with  u(z)  the  total  acoustic 
field,  that  is: 


|  The  research  reported  in  this  paper  has  been  made  possible  through  the  support  and  sponsorship  of  the  Italian 
Government  through  the  Mimstero  per  PUniversita  e  per  l&  Ricerca  Scientific*  under  contract  MURST  40% 
1990  and  of  the  U.  S.  Government  through  the  Air  Force  Office  of  Scientific  Research  under  contract  n. 
AFOSR  90*0228 
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u(i)  =  u'(x)  +  u’(i)  (1.2) 

The  total  acoustic  field  u(r)  satisfies  the  Helmholtz  equation: 

Au  +  fc2u  =  0  in  R3\D  (1.3) 


and  the  scattered  acoustic  field  u*(x)  satisfies  the  Sommerfeld  radiation  condition  at 
infinity,  that  is: 


„  Ijm  i|x[| 

!|rtf  — oo 


.  du 3 


-  iku  M  =  0 


(1-4) 


where  A  =  -j-;  +  j b.  -t-  is  the  Laplace  operator.  Moreover  the  total  acoustic  field 
u(r)  satisfies  a  boundary  condition  on  dD.  This  boundary  condition  can  be  formulated 
in  several  different  ways,  depending  on  the  nature  of  the  obstacle  D . 

In  [2], [3], [4]  we  have  considered  the  acoustically  soft  obstacles  that  are  characterized  by 
the  Dirichlet  boundary  condition: 

u  =  0  on  dD  (1.5) 


In  [4],(5j  we  have  considered  the  acoustically  hard  obstacles  characterized  by  the 
Neumann  boundary  condition 


—  =  0  on  ac 
au 


(1.6) 


where  i>  is  the  unit  normal  on  dD,  and  the  obstacles  characterized  by  an  acoustic 
impedance  that  satisfy  the  mixed  boundary  condition 

u  +  X^--=0  on  dD  (1.7) 

c w 

We  assume  that  x  is  a  given  constant.  We  consider  three  boundary  value  problems:  the 
Dirichlet  boundary  value  problem  given  by  (1.3),(1.4),(  1.5),  the  Neumann  boundary 
value  problem  gi  'cn  by  (1.3),  (1.4), (1.6)  and  the  mixed  boundary  value  problem  given 
by  (1.3),  (1.4), (1.7).  In  (6J  it  is  shown  that  the  scattered  field  ti  ’(x)  of  the  Dirichlet,  the 
Neumann  and  mixed  boundary  value  problem  has  the  following  expansion 


<*’(*)  =  -  rF0(z,k,a)  +  0(  — 

„A|  '1*11 


■ )  when  ||ij[  — •  oo 


(1-8) 


where  £  =  ,  x  0  and  F0(£,h,a)  is  the  far  field  pattern  generated  by  the  incoming 

wave  (1.1)  that  hits  the  obstacle  D. 

In  (2I,(3|,[4|,[5|  we  have  introduced  a  numerical  method  for  an  inverse  problem  for  the 
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three  dimensional  Helmholtz  equation,  that  is  from  the  knowledge  of  the  nature  of  the 
obstacle  i.e.  the  boundary  condition  on  dD  satisfied  by  u(z)  and  from  the  far  fields  F0 
generated  by  several  incoming  waves  we  want  to  recover  the  shape  of  the  obstacle  dD. 
To  be  more  precise  let  An  ,  n  =  1,2,,..  be  the  eigenvalues  of  the  Laplace  operator  in 
the  interior  of  D,  with  Dirichlet  boundary  condition  (1.5)  or  with  Neumann  boundary 
condition  (1.6)  or  with  the  mixed  boundary  condition  (1.7);  let  B  =  (r  £  R3  |  ||zj|  <  1}, 
and  dB  be  the  boundary  of  B.  We  will  consider  the  following  inverse  problem: 

Problem  1.1  /ntierje  acoustic  probi'-'n.  Let  us  assume  that  u(z)  satisfies  the  Dirichlet 
boundary  condition  (1.5)  or  the  Neumann  boundary  condition  (1.6)  or  the  mixed 
boundary  condition  (1.7).  Let  fl]  C  dB  ,  fii  C  {(  £  R  |  f  >  0}  be  two  given  sets 
such  that  A,  £  Qi  i  —  1,2...  From  the  knowledge  of  F0[x,  k,  a),  for  a  6  fli  ,  —  k2  g 
determine  the  boundary  of  the  obstacle  dD. 

We  note  that  the  condition  A,  g  1  —  1,2...  is  a  non-resonance  condition,  that  Qj  is 
the  set  of  the  directions  of  the  incoming  waves  and  that  the  far  fieid  F0  is  observed  for 
x  g  dB.  Let  us  now  consider  the  electromagnetic  problem.  For  time  harmonic  waves 
the  Maxwell  equations  are  reduced  to  the  time  harmonic  Maxwell  equations  [6]  (see 
chapter  4). 

In  the  following  we  wiil  use  occasionally  complex  vectors  abusing  of  the  notations. 

Let  E'(x)  be  the  electric  field  associated  to  a  linearly  polarized  time  harmonic  incoming 
wave  propagating  in  an  homogeneous  isotropic  medium,  that  is: 

E'(x)  =  we,‘(-l'2)  (1.9) 

where  to,  a  €  R3  with  ||a||  =  1  are  given  and  k  >  0  is  the  wave  number,  moreover  we 
assume  that: 

divE'(x)  =  ik(w,  a)e'*f--)  =  0  (1.10) 

where  E’(x)  =  (£;(z),  £j(x),  £l(z))  and  divE'(x)  =  ^  We 

note  that  to  is  the  polarization  vector  and  a  is  the  propagation  direction  of  the  incomig 
electric  field.  We  note  that  the  magnetic  field  H'(x)  associated  to  this  incoming  wave 
is  given  by: 

&‘{x)  =  IrCurlE'ix)  (1.11) 

ik 


where  curLE'(z)= 

Let  us  denote  with  E’ (x) 


_  3 £',(£)  dEUl)  _  3E‘,(r)  U)  _  3g;(rA 

dz  ’  dz  dt  ’  dz  3jf  j  ' 

the  electric  field  scattered  by  the  obstacle  D  when  hitted  by 


the  incoming  wave  E'< x)  and  with 


E(x)  =  E\x)  +  E’{x) 


(1.12) 


the  total  electric  field.  It  is  easy  -o  see  ,’6J  (see  chapter  4)  that  the  time  harmonic  Maxwell 
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equations  in  an  homogeneous  isotropic  medium  that  does  not  contain  eiectric  charge 
reduce  to  the  vector  Helmholtz  equation,  so  that  the  scattered  field  E_’(z)  satisfies: 

&E’(x)  +  k2E:'{x)  =  0  in  R3  \  Z3  (1.13) 


together  with  the  equation: 

divjd'is.)  =0  in  R3  \  D  (1.14) 

while  the  Silver-Muiler  radiation  condition  at  infinity  [6]  (see  paragraph  4.2)  reduces  to: 


curlE’(x)  xi  -  ikE'(x)  =  o 


ll£il  00 


(1.15) 


where  x  is  the  vector  product.  Let  v(x)  be  the  exterior  unit  normal  to  dD ,  for  a 
perfectly  conducting  obstacle  D,  we  will  assume  the  following  boundary  condition: 

£(x)  x  kis.)  —  Q  .  x  €  dD  (1-16) 


In  [6]  it  is  shown  that  E'{x ),  solution  of  the  boundary  value  problem  (1.13),  (1.14), 
(1.15),  (1.16)  has  the  following  expansion: 

e't|l*ll  /  i  \ 

£’U)  =  -jj^ip^o (i,  k,a,u>)  +  0  J  ,  IUH  -  oo  (1.17) 

where  E0{x,k,a,w)  is  the  (electric)  far  field  pattern  generated  by  the  incoming  wave 
(1.9)  that  hits  the  obstacle  D. 

In  (7j  we  introduce  a  numerical  method  for  an  inverse  problem  for  the  three  dimensional 
vector  Helmholtz  equation  (1.13).  That  is,  from  the  knowledge  of  the  nature  of  the 
obstacle  (i.e.  the  fact  that  the  obstacle  is  perfectly  conducting)  and  from  the  (electric) 
far  fields  E0(x,  k,  <2,  tya)  generated  by  several  (known)  incoming  waves  we  want  to 
recover  the  shape  of  the  obstacle  dD. 

To  be  more  precise  let  A„  n  =  1,2, ...  be  the  eigenvalues  of  the  "vector”  Laplace 
operator  restricted  to  the  divergence  free  vector  fields  with  the  homogeneous  boundary 
condition  (1.16)  in  the  interior  of  D.  In  [7]  we  have  considered  the  following  inverse 
problem: 

Problem  1.2  Inverse  electromagnetic  problem.  Let  f2t  C  dB  ,  Qj  C  R3  ,  fij  C  R  be 
three  given  sets  such  that  \,  £  f23  t  =  1,2, ...  and  let  E0(x,k,a,wa)  be  the  (electric) 
far  field  defined  in  (1.171.  From  the  knowledge  of  £0(r,  k,  a,  wa)  for  agfi ,,  wa  G  ftj, 
—  k2  €  H}  determine  the  boundary  of  the  obstacle  dD. 

We  note  that  the  condition  \,  i  '.1,  i  =  1,2, . . .  is  a  non  resonance  condition,  that  f2|  is 
the  set  of  the  directions  of  the  incoming  waves  and  that  122  is  the  set  of  the  polarization 
vectors  of  the  incomig  wav*>s 
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The  inverse  Problems  1.1,  '..2  have  been  studied  generalizing  the  "Herglotz  function 
method"  introduced  by  CT  ton  and  Monk  [1|.  This  technique  is  supposed  to  be 
particularly  effective  in  the  resonance  region,  that  is,  when 

IcL  ~  1  '  (1.18) 

where  L  is  a  characteristic  length  of  the  obstacle  D. 

In  section  2  we  develop  the  mathematical  relations  needed  to  develop  our  procedure  and 
we  outline  the  numerical  method  derived  from  them.  Finally  in  section  3  we  present 
some  numerical  results. 

2.  THE  MATHEMATICAL  FORMULATION  OF  THE  INVERSE 
PROBLEM 


We  will  restrict  our  attention  here  to  the  inverse  Problem  1.1  associated  to  the 
Dirichlet  boundary  value  problem  (1.3), (1.4), (1.5).  The  remaining  inverse  problems  are 
treated  in  a  similar  way  and  we  refer  to  (3j,(4[,(5j,(?j. 

For  x,y  €  R3  let 


*(*lla-2ll) 


Ws.  -  yll 


(2.1) 


be  the  Green's  function  of  the  Helmholtz  operator  with  the  Sommerfe’.d  radiation 
condition  at  infinity.  It  is  easy  to  see  that: 


<&(fc||x-y||)  =  ^-e-,*(1'2)+0(^)  when  ||*||-oo 

moreover  from  the  Helmholtz  formula  [6]  we  h.ave: 

1  ft  t  X  if  xeD 

JaD  ~  <My)  <My)  a'  ~  l  «*(*)  if  x  e  R3  \  d 


where  dcri(y)  is  the  surface  measure  on  dD. 
Substituting  (2.2)  in  (2.3)  and  using  (1.8)  we  have: 


F0{x,  k,a) 


_1_ 
4  it 


5e-.k(£.y) 

du(y) 


3«(y) 

a„(y)e 


■*(i'.Jf))<f<r,(y) 


(2.2) 


(2.3) 


(2.4) 


Let  g(x)  €  L2(3B,di72)  where  L2(3B,d<72)  is  the  space  of  square  integrable  functions 
with  respect  to  the  surface  measure  on  the  unit  sphere  d<Ji  and  g{i )  is  the  complex 
conjugate  of  y(x). 

For  every  g(x)  €  L7(dB,  do-i)  from  (2.4)  interchanging  the  integrals  we  have: 
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I  F0{i,K,a):j(i)da2{i)=  I 
JdB  -lx  Jgg 

f  ,  ,  du(y)  k(i 


Vao  -  a  V(y)  dvyy) 

1  f  .  ,  dv(ky)  du{j), 

=  7-  /  "1/  777  _  7-7=7 v{ky))d<7X(y) 

4 *  JaD  -  3i/(y)  oi'(y) 


where 


*(*y)  =  /  y(x)e'*(i'ii,^2(x) 

JdB 


It  is  easy  to  sec  that  v(ky)  satisfies  the  Helmholtz  equation  for  y  g  R3,  moreover  v(ky') 
is  the  Herglotz  wave  function  corresponding  to  the  Herglotz  kernel  y(i).  Since  the  total 
acoustic  field  u  satisfies  the  boundary  condition  (1.5)  on  the  surface  dD  formula  (2.5) 
reduces  to: 


/  F*(k,  k,a)g(i)do2{x)  =  f 
JdB  4T  Jg 


1  f  du(y) - 

—  /  y-p-v(fcy)ci<7,(y) 

Jao  dt/{y) 


We  restrict  our  attention  to  the  Dirichlet  Herglotz  domains,  that  is  domains  such  that 
the  unique  solution  v  of: 

(A  +  &2)v  =  0  x  €  D  (2.8) 


V 


"fciisr 


x  g  OD 


(2.9) 


is  given  by  (2.6)  for  a  suitable  choice  gn(i.}k)  of  g(£) 

We  note  that  in  the  definition  of  Herglotz  domains  we  have  exploited  the  hypothesis 
-Jfc2  A,  i  =  1,2,  •••  .  A  simple  computation  shows  that  the  sphere  with  center  the 
origin  is  an  Herglotz  domain  that  is  the  class  of  the  Herglotz  domains  is  not  empty.  In 
(2.7),  let  v  be  the  Herglotz  wave  function  associated  to  gn(i,k ),  using  (2.9)  and  the 
Helmholtz  formula  (2.3)  we  have: 


S.,k,a)gu(i,k)drj2(i)  =  - 
k 


(2.10) 


when  gu(£,k)  is  the  Herglotz  kernel,  formula  (2.10)  holds  'ik.a.  Problem  1.1  proposed 
in  section  1  will  be  solved  in  three  steps: 

(i)  from  the  knowledge  of  some  far  fields  F0  using  (2.10)  determine  the  Herglotz  kernel 
gn(x,k)  of  the  domain  D 

(ii)  from  gn{x,k)  using  (2.6)  find  the  corresponding  Herglotz  wave  function  i; 
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(iii)  determine  dD  using  (2.9) 

The  steps  (i),(ii),(iii)  are  performed  with  a  fixed  k,  however  when  necessary  several  \  two 
or  three)  values  of  k  can  be  used  [4),[5|. 

The  numerical  method  that  exploits  the  previously  described  reconstruction  procedure 
consists  of  some  linear  algebra  computations  to  perform  step  (i),  explicit  analytic 
computation  to  perform  step  (ii),  a  global  minimization  procedure  to  perform  step 
(iii).  The  global  minimization  algorithm  used  is  based  on  the  numerical  integration  of 
a  system  of  stochastic  differential  equations  [9|,(10).  Further  details  are  contained  in 
[2U3],[4),(5U7|. 

3.  NUMERICAL  RESULTS 


The  surfaces  dD  considered  are  the  following  ones: 


1- 

Oblate  Ellipsoid 

(f1)2  +  (§y)2  4-  z2  =  1 

(3.1) 

2- 

Prolate  Ellipsoid 

x2  +  y2  +(|z)2  =  1 

(3.2) 

3- 

Vogel's  Peanut 

r  =  |(cos2  0  +  |  sin2  8)]/2 

(3.3) 

4- 

Short  Cylinder 

((§x)2+(fy)2)5  +  z>°  =  l 

(3.4) 

5- 

Pseudo  Apollo 

r  =  |(!^  +  2  cos  30)' f2 

(3.5) 

All  these  surfaces  are  cylindrically  symmetric  with  respect  to  the  z-axis  and  the  surfaces 
1,2, 3, 4  are  also  symmetric  with  respect  to  the  equator. 

We  observe  that  the  obstacles  D  corresponding  to  1,2,4  are  convex  and  the  ones 
corresponding  to  3,5  are  not  convex.  Finally  a  characteristic  length  L  of  the  obstacles 
can  be  chosen  equal  1. 

Let  ( )  be  polar  cohordinates  and  i  =  (sin 0 cos  sin# sin  <£, cost?),  dD  =  (r  = 
/(0,  0)|O<0<t,O<^<  2zr }  for  j  =  0, 1, ...  ,36  let  9,  =  f(9j,  0)  be  the  exact 

values  of  the  surfaces  given  by  (3.1),. . -,(3.5)  and  /c(0;,O)  be  the  values  reconstructed 
performing  the  numerical  procedure  described  in  section  2-  The  relative  L 2  error  in  the 
points  {(0j,O)  |  j  =0,1,...  ,36},  that  is  : 

£(/(0;.O)-/c(0„O))2|  3 

El,  =  - - -  (3.6) 

Em.o) 

l  J=° 

is  used  as  a  performance  index. 

In  the  results  shown  below  we  use  9-11  different  directions  for  the  incoming  waves  and 
in  the  case  of  the  electromagnetic  problem  two  linearly  independent  polarizations  for 
each  incoming  direction.  The  far  fields  are  obtained  from  the  numerical  solution  of 
the  corresponding  boundary  value  problems.  That  is,  for  example,  for  the  acoustic 
Dirichlet  inverse  problem  the  boundary  value  problem  ( 1.3),(  1.4),(  1.5).  A  random  error 
of  order  l%-5%  is  added  to  the  numerically  computed  far  fields  in  order  to  simulate 
actual  measurements.  The  far  fields  are  supposed  to  be  known  on  the  full  solid  angle 
(complete  data)  or  only  on  a  finite  set  (9-11)  directions  (incomplete  data). 


47 


The  resuits  obtained  with  our  reconstruction  procedure  are  •'  '»n  in  tables  3.1,  3. 
3.3,  3.4. 


Table  3.1 

Acoustic  Dirichlet  inverse  problem 


Object  |  k 

.N umber  of 
incoming  wivfj 

Complete 

data 

El* 

Oblate  Ellipsoid!  3 

9 

no 

0.0071 

Prolate  Ellipsoid!  3 

9 

no 

0.031 

Vogel’s  Peanut  |  3  j  9 

yes 

0.019 

Short  Cylinder  |  3 

9 

no 

0.0091! 

Pseudo  Apoilo  (  4 

9 

yes 

0.062 

Table  3.2 


Acoustic  Neumann  inverse  problem 


Object 

LL 

Number  of 
incoming  waves 

Complete 

data 

Oblate  Ellipsoid 

3 

9 

no 

0.0079 

Prolate  Ellipsoid 

3 

9 

no 

0.0139 

Vogel’s  Peanut 

3 

9 

no 

0.0463 

Short  Cylinder 

2 

9 

yes 

0.0476 

Pseudo  Apollo 

2-3 

9 

yes 

0.0353 

Table  3.3 


Acoustic  Mixed  inverse  problem  (x  =  1) 


Object 

D 

Oblate  Ellipsoid 

D 

9 

yes 

Prolate  Ellipsoid 

4 

9 

yes 

0.017 

Vogel’s  Peanut 

El 

9 

no 

imwrt 

Short  Cylinder 

11 

9 

yes 

0.014 

Pseudo  Apollo 

2-3 

9 

yes 

0.036 

Table  3.4 


Electromagnetic  inverse  problem 


Object 

k 

.Number  ol 
incoming  waves 

Complete 

data 

El> 

Oblate  Ellipsoid 

3 

is 

yes 

fflKETi? 

Prolate  Ellipsoid 

3 

18 

yes 

0.0323 

Vogel’s  Peanut 

3 

IS 

riTiVI&i 

ISAiCaik 

Short  Cylinder 

2.5 

22 

MtKttl 

Pseudo  Apoilo 

22 

yes 

0.0411 

As  an  example  in  the  figures  1  and  2  the  quality  of  our  reconstructions  is  shown. 
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ORIGINAL  RECONSTRUCTED 


Short  Cylinder 
ORIGINAL 


Pseudo  Apollo 
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